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ABSTRACT 

We present an update of our low-resolution spectroscopic follow-up and model atmosphere 
analysis of hot subdwarf stars from the Galaxy Evolution Explorer (GALEX) survey. Targets 
were selected on the basis of colour indices calculated from the GALEX GR6 N\j\ , Guide Star 
Catalogue (GSC2.3.2) V and the Two Micron All Sky Survey (2MASS) J and H photometry. 
High signal-to-noise ratio spectra were obtained at the European Southern Observatory (ESO) 
and the Kitt Peak National Observatory (KPNO) over the course of three years. Detailed H, 
He and CNO abundance analysis helped us improve our r e ff , log g and He abundance deter- 
mination and to constrain CNO abundances. We processed 191 observations of 180 targets 
and found 124 sdB and 42 sdO stars in this sample while some blue horizontal branch stars 
were also found in this programme. With quantitative binary decomposition of 29 composite 
spectra we investigated the incidence of A, F and G type companions. The incidence of late 
G and K type companions and their effects on subdwarf atmospheric parameters were also 
examined. 

Key words: catalogues - surveys - stars: abundances - stars: atmospheres - binaries: spec- 
troscopic - subdwarfs 



1 INTRODUCTION 

Hot subdwarfs are core He burning stars located at the blue end 
of the horizontal branch (HB), also known as the extreme horizon- 
tal branch (EHB). These subluminous objects are located roughly 
midway between the early-type main-sequence (MS) stars and the 
white dwarf (WD) cooling curve in the Hertzsprung-Russell dia- 
gram (HRD). Being in a relatively long lasting (~ 160 Myr), inter- 
mediate evolutionary stage of ~1M Q stars, they are quite common 
and likely to overwhelm WDs in blue and ultraviolet (UV) surveys 
of old stellar populations. These stars are the primary sources of 
the UV excess of elliptical galaxies (Brown et al. 1997, O'Connell 
1999). 

To look for bright, thus nearby WD and subdwarf candidates 
in the Galaxy Evolution Explorer (GALEX) data base, Vennes, 
Kawka & Nemeth (2011a) (hereafter: Paper I) devised a method 
based on UV, optical and infrared (IR) colours. Between 2008 
and 2011, our systematic, low-resolution spectroscopic follow- 
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up with the European Southern Observatory (ESO)/New Tech- 
nology Telescope (NTT) and the Kitt Peak National Observatory 
(KPNO)/Mayall telescopes confirmed 166 stars as hot subdwarfs. 
Here, we present the results of a model atmosphere analysis of 180 
stars. Details of previous surveys of blue or ultraviolet excess ob- 
jects and of our source selection can be found in Paper I. Here we 
modified the original selection criteria to avoid de-selecting bright 
sources with unreliable GALEX Fu\ photometry. For GALEX 
N\jy < 12.5, where iVuv is corrected for non-linearity effects (Mor- 
rissey et al. 2007), we selected targets with Nu\ — V < 0.5, where 
V is the Guide Star Catalogue version 2.3.3 (GSC2.3.2; Lasker et 
al. 2008) photographic magnitude (see Paper I), while ignoring the 
GALEX Fuv magnitude. 13 bright objects were added to our sam- 
ple by using the modified criteria. Among the 1 80 stars we found 
124 subdwarf B (sdB) and 42 subdwarf O (sdO) stars. Most of the 
remaining objects are classified as blue HB (BHB) stars. The 52 
stars from Paper I are also re-analysed here with an extended abun- 
dance pattern and using a new fitting algorithm that provided con- 
sistent results with Paper I. Such large and homogeneously mod- 
elled samples of bright subdwarfs are useful in identifying candi- 
dates for pulsation (e.g., Ostensen et al. 2010), radial velocity stud- 
ies (e.g., Geier et al. 201 la, Napiwotzki 2008) and eventually to set 
observational constraints for formation and evolution theories. 

The present paper deals with optical/IR data exclusively. We 
list A'uv magnitudes, but refer the interested reader to the UV work 



2 P. Nemeth, A. Kawka and S. Vermes 



of Girven et al. (2012). A similar UV analysis will be reported else- 
where. 

There are two main formation concepts for subdwarf stars, but 
current observational data cannot satisfyingly distinguish the con- 
tribution of their formation channels. The canonical subdwarf for- 
mation theory (Mengel, Norris & Gross 1976) invokes binarity and 
strong mass loss on the red-giant branch (RGB) for the formation 
of sdB stars, and WD binary mergers for the formation of He-rich 
hot sdO (He-sdO) stars (Webbink 1984). For a recent review of the 
canonical formation channels and binary population synthesis see 
Hu et al. (2008), Han et al. (2003a) and Han et al. (2003b). 

The late hot-flasher scenario (D'Cruz et al. 1996) suggests a 
different evolutionary path that explains the abundance diversities 
observed in hot subdwarfs. This model also assumes a strong mass- 
loss on the RGB, but the star undergoes a core He-flash on the WD 
cooling track. Then, convective mixing would bring core material 
into the hydrogen envelope and would account for He, C and N 
overabundances. 

Both models require an effective process to remove most of 
the hydrogen envelope on the RGB. The main difference is that 
the canonical formation scenario assumes a core flash on the tip of 
the RGB without mixing, while in the late hot-flasher scenario the 
star departs from the RGB before the flash happens allowing for 
shallow or deep convective mixing (Lanz et al. 2004). 

Stringent observational constraints are needed to determine 
the relative contribution of these formation models. For this reason, 
detailed model atmosphere and binary analyses of large samples 
carried-out independently with different model atmosphere codes 
are necessary to reveal the contribution of various formation chan- 
nels and the condition of mass loss on the RGB. 

Throughout the paper we use the notation [He/H]. = 
log l0 (« H e/' 1 H) f° r abundance and [Fe/H] = [Fe/H]. - [Fe/H] Q for 
metallicity. We omit the GALEX prefix, as well as the seconds and 
arcseconds of the identifier when referring to stars in the catalogue. 
Finally, we refer to possible short- (p-mode) and long-period (g- 
mode) pulsators as rapid and slow pulsators using the nomenclature 
of Kilkenny et al. (2010). 

In Section 2 we briefly describe our data acquisition and re- 
duction. In Section 3 we give details of our spectral modelling and 
fitting with TLUSTY and XTGRID, and our approach for binary de- 
composition. In Section 4 we present the properties of hot subd- 
warfs in the T eff - log g and r cff - He diagrams, and describe our 
results on CNO abundance trends and binarity. In Section 5 we in- 
vestigate the sample completeness by comparing our observed lu- 
minosity distribution to a modelled one. Modelling homogeneity is 
examined in Section 6. Subdwarf populations are reviewed in light 
of our results in Section 7. The catalogue of our GALEX sample 
is presented in Section 8 and Section 9 summarizes our results. 
Finally, in Section 10, we propose to initiate a large collaborative 
study on hot subdwarf stars. 



2 FOLLOW-UP SPECTROSCOPY 

Low-dispersion optical spectra were obtained at two sites during 
seven observing runs between 2008 and 201 1. At ESO, we used the 
ESO Faint object Spectrograph and Camera (EFOSC2) attached to 
the 3.6-m NTT at the La Silla Observatory. On UT 2008 October 
19-22, 2009 March 2-4 and 2009 August 23-27, we used grism#ll 
(300 lines mm" 1 ) with a dispersion of ~4.17 A per binned pixel (2 
x 2). With a 1 arcsec slit width we obtained a resolution of AA « 
13.7 A. On UT 2010 March 2-4 and 2010 September 18-21 we 



employed grism #7 (600 lines mm ~ 1 .96 A per binned pixel) 
and a 1 arcsec slit width resulting in a resolution of AA « 6.4 A. 

At the KPNO on UT 2010 March 23-26 and 201 1 January 28- 
31 the Ritchey-Chretien Focus Spectrograph (RC Spectrograph) 
was used with the 4-m Mayall telescope and KPC-10A (316 lines 
mm" 1 ) grating, delivering a spectral resolution of AA « 5 A in first 
order. The slit width was set at 1.5 arcsec resulting in a dispersion 
of ~2.75 A per pixel. 

The KPNO spectra and the lower-resolution ESO data cover 
the entire Balmer-series from 3700 to 7400 A. The ESO AA a 6.4 
A data covers the upper B aimer line series from 3600 to 5200 A. 
All data were reduced using standard IRAF 1 procedures. 

Our selection also included three hot subdwarf candidates pro- 
posed by Jimenez-Esteban, Caballero, & Solano (2011). Two of 
these objects are GALEX sources outside of our selection crite- 
ria (TYC 6017-419-1 and TYC 9327-1311-1), while another lies 
in a region not covered by the survey (TYC 9044-1653-1). We in- 
cluded another two targets (J0716+2319 and J2349+41 19) that are 
not subdwarfs, but accidentally met the criteria due to their erro- 
neous V magnitude in our source selection. Although these stars 
are not formally part of the UV selection, we analyze their spectra 
in Section 3. 



3 SPECTRAL ANALYSIS 

3.1 Model atmospheres 

We computed H/He/CNO non-LTE model atmospheres with 
TLUSTY 200 and synthetic spectra with SYNSPEC 48 (Hubeny & 
Lanz 1995; Lanz & Hubeny 1995). Model atoms of Hi, Hei-n, 
Cn-iv, Nm-v and Oiv-vi; and detailed line profiles of H and He 
were used in SYNSPEC. For all relevant ions, we included the most 
detailed model atoms from the OSTAR2002 (Lanz & Hubeny 2003) 
and BSTAR2006 (Lanz & Hubeny 2007) data base. 

Model atmospheres for BHB stars with effective temperatures 
below 20 000 K were calculated with H i, He i-n, C n-iv, N ii-iv and 

n-m model atoms. A line list compiled from Kurucz CD-ROM 
23 and available at the SYNSPEC web page 2 was used. Details of 
the model atoms used in this work are listed in Table 1 . 

3.2 Spectral fitting with XTGRID 

A spectral fitting method was employed by combining the steepest- 
descent and simplex algorithms, implemented in our new x 1 min- 
imizing fitting program XTGRID. This Python program is an ad- 
justable interface for TLUSTY and SYNSPEC, designed to carry 
out iterative multi-wavelength spectral analysis of hot stars from 
soft X-rays to near infrared wavelengths. The procedure requires 
a starting model and with successive adjustments approaches the 
observed spectrum. 

At the beginning, XTGRID follows the gradient of steepest de- 
scent with a maximal step size that TLUSTY convergence allows. 
New step directions and sizes are calculated in a small grid (sim- 
plex step) and adjusted independently for each parameter. To follow 
the steepest gradient, relative step sizes are normalized to the lead- 
ing parameter. Model atmospheres and synthetic spectra are calcu- 
lated for a new temperature and gravity in all iterations, and only 
synthetic spectra are calculated for new abundances. Directions are 

1 http://iraf.noao.edu/ 

2 http://nova.astro.umd.edu/ 
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Table 1. Model atoms used for hot (H) and cool (C) subdwarfs. 



Ion Number of levels Frequency of Model 
and line transitions highest level, s 
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calculated from these fits and the starting model gets updated with 
the new model. This process is pursued until all relative changes 
(T, logg, abundances and^ 2 ) decrease below 0.5 per cent in three 
consecutive iterations and all maximum allowed step sizes decrease 
below 50 per cent of their initial values. The latter condition ensures 
that a minimum is found for parameters with small gradients. 

Temperature and gravity have a higher convergence rate than 
abundances. To take advantage of this observation and accelerate 
our procedure, parameter relaxation is included in the program. 
When relative changes of temperature or gravity decrease below 
0.5 per cent they can be kept fixed for five iterations. This allows 
for slow converging parameters (like abundances) to approach their 
final values faster. Furthermore, such iterations save two out of the 
three model atmosphere calculations per iteration, the most time 
consuming part. After 1 1 iterations all models are calculated again 
for three iterations, to examine relative changes and to test if the 
final model is achieved. If the convergence limit is not reached, 
this loop of 14 iterations is repeated. This strategy provides about a 
factor of 2 speed-up with respect to regular iterations when model 
atmospheres are always calculated for a new temperature and grav- 
ity. Additionally, the fitting procedure accelerates automatically as 
by approaching the solution the relative differences between mod- 
els get smaller. 

XTGRID regularly reads the optional TLUSTY flag file, in 
which important control parameters can be stored for the fitting 
process. This way it is possible to adjust any fitting parameters or 
procedures during calculations. 

A critical point in spectral fitting is the method how models are 
normalized to observations. Over a wide spectral range the match 
between theoretical and observed continua is usually poor in the be- 
ginning, while over a short wavelength range the continuum might 
not be visible because of broad spectral lines or line blanketing. 
To account for both effects, XTGRID samples the spectrum in small 
sections. The selection of these fitting windows is either specified 
by the user or the program partitions the entire data in equally sized 



sections. The first setup can be used for selecting certain lines or 
ranges for fitting, like the hydrogen Balmer series or certain metal 
lines. When the whole spectrum is considered, only the segment 
size needs to be specified. The optimal section size depends on the 
spectral resolution, the spectral features being analysed and must 
be selected carefully. The synthetic spectrum is then normalized 
in each of these bins to the observed data by finding the ratio of 
their median fluxes. This simple approach tends to normalize for 
the continuum and normalization improves with the goodness of fit. 
Low order variations of the spectrum are significantly reduced this 
way, therefore the effects of interstellar reddening and poor fluxing, 
often seen at short optical wavelengths, are treated. An additional 
benefit of this method is, when the whole spectrum is used, that the 
X 1 considers the interconnections between all spectral lines and the 
continuum. This can be very helpful in investigating the effects of 
metal abundances on hydrogen lines, in particular, to address the 
Balmer-line problem. 

Although low order variations can be significantly reduced 
with a careful sectioning of the spectrum, our normalisation method 
considers the continuum as well, hence it is sensitive for flux cal- 
ibration. Large deviations from the theoretical continuum due to a 
poor flux calibration can affect the final parameters. Hence, for con- 
sistent modelling the same sampling of the spectrum is necessary, 
or to work with continuum normalized spectra. XTGRID is designed 
to fit continuum normalized data as well. 

To match the resolution of synthetic spectra to the observa- 
tions, they are convolved with a Gaussian profile at either con- 
stant resolution (full width at half-maximum) or constant resolving 
power R = A/AA. 

Radial velocity corrections can be defined by specifying AA/A 
or taken as free parameter. If fitted, a^ 2 is calculated at small shifts 
of the model and the minimum is determined by a low order poly- 
nomial fit. This correction yields absolute radial velocities and up- 
dated after every third iteration with adjusted shifts. 

After fitting is done, parameter errors for 60, 90 and 99 per 
cent confidence intervals are estimated by mapping the Ax 2 with 
respect to the final abundance (X) at representative points in the 
range of -0.99 < AX/X < 100. Temperature and gravity errors 
are measured similarly in the range of -0.26 < AT IT < 0.26 and 
-0.195 < A logg/ logg < 0.195. Parameters are changed until the 
statistical limit for the 60 per cent confidence level at the given 
number of free parameters is reached. Model atmospheres are cal- 
culated during error calculations for all synthetic spectra to track 
small changes of the^ 2 . This way the size and symmetry of error 
intervals can reflect the quality of the parameter determination. Er- 
rors at 90 and 99 per cent are extrapolated by a parabolic fit for the 
upper and lower error intervals independently. 

XTGRID is scalable for cluster calculation and to further accel- 
erate the fitting procedure, previously calculated models are reused 
until relative changes of temperature and gravity drop below 13 per 
cent of their respective maximum allowed values. Following this 
procedure, the program builds a model cloud that is more closely 
spaced near best solutions. Such accelerations are necessary to cope 
with the computing demand of model atmosphere analyses. In par- 
ticular, in batch analysis of multiple spectra, because any change in 
the input data, physics or fitting procedure requires the recalcula- 
tion of all fits to maintain homogeneity. In the final step LTpX fit 
summaries and Gnuplot scripts are produced for graphs shown in 
this paper. Figure 1 shows the best-fitting of the He-rich BHB star 
J2059-4232. 
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Figure 1. Observed NTT and best-fitting model spectrum for J2059^4232. Identifications are given for He and CNO lines with theoretical equivalent widths 
larger than 1 A. The presence of He, C and N is statistically confirmed. 



3.3 Binary spectral decomposition 

Binarity plays an important role in the theory of subdwarf for- 
mation and evolution. In the canonical formation theory (Han et 
al. 2003a) the common-envelope, Roche lobe overflow and binary 
merger channels all depend on binary parameters. About 40 per 
cent of sdBs were found in binary or multiple systems with MS 
companions (Maxted et al. 2001, Napiwotzki, Karl & Lisker 2004). 
Reed & Stiening (2004) derived an sdB-MS (AO to M2) binary 
fraction of 53 ± 6 per cent by using Two Micron All Sky Survey 
(2MASS; Skrutskie et al. 2006) J - H and optical B - V colours. 
Thejll, Ulla & MacDonald (1995) found ~43 per cent composite 
spectra in their sample of 27 subdwarfs. Stark & Wade (2003) 
found 40 per cent composite spectra in their magnitude limited 
sample of about 600 subdwarfs collected from the Palomar-Green 
survey (PG, Green, Schmidt & Liebert 1986) and the Catalogue 
of Spectroscopically Identified Hot Subdwarfs (Kilkenny, Heber & 
Drilling 1988). At least 20 per cent of subdwarfs show significant 
IR excess in the V—J - J-H diagram (Ferguson, Green & Liebert 
1984). In agreement with this number, about 19 per cent of our sam- 
ple has V—J > and J-H > indicative of an IR excess (Figure 
2). The true fraction of composites depends heavily on the quality 
of the spectra [signal-to-noise ratio (SNR), resolution] and selec- 
tion criteria. Combining the frequency of composite spectra, radial 
velocity variables and photometric variables (ellipsoidal variations) 
we conclude that the binary fraction is close to two-thirds in agree- 
ment with Maxted et al. (2001). 

About 17 per cent of our sample shows strong double lined 
composite spectra and about 33 per cent have noticeable Can H&K 
013933 and ,13968 A), Mgi 015183 A) or Nai D (,15890 and A5896 
A) absorption lines. Fe i and CH molecular lines blend at ~/14300 
A and notable in F and G type stars. This strong feature is often re- 
ferred in the literature as the G band. We found G-band absorption 
in composite spectra. Although, these lines might have an interstel- 
lar origin, their simultaneous presence and a flat flux distribution 
can indicate a composite spectrum. 

The high incidence rate of binaries and their importance in 
subdwarf evolution in addition to the fact that subdwarfs in com- 
posite spectra binaries cannot be modelled as single stars, inspired 



us to include binary decomposition in XTGRID. In double-lined bi- 
naries both components can be examined simultaneously. In partic- 
ular, the F and G companions are relatively easy to identify because 
they have characteristic spectral features different from subdwarfs 
and a comparable optical brightness. Late G and K companions 
present a bigger challenge because they have significantly lower 
contribution and weaker lines. 

A library of low- and medium-resolution observed spectra was 
collected for spectral decomposition from the MILES (Cenarro et al. 
2007, Sanchez-Blazquez et al. 2006) and HILIB (Pickles 1998) li- 
braries with resolutions AA = 2.3 A and R = 500, respectively. For 
visual inspections the spectral library of Silva & Cornell (1992) was 
also used. Altogether, 946 spectra between ,13525 and ,17500 A at 
2.3 A resolution were included in our decomposition library from 
the MILES data base. This library provides a wide range of spectral 
types at sufficiently high resolution and good fluxing. HILIB tem- 
plates were used to calculate synthetic colour indices in the optical 
and infrared. 

We searched for the best-fitting template spectrum by interpo- 
lating in temperature, surface gravity, metallicity (as [Fe/H]) and 
flux ratio (F, )iCO „, p J(F A comp + F A sd )) of the secondary component 
along with TLUSTY model parameters for the primary. The decom- 
position was performed at every third iteration and during the three 
consecutive iterations when fitting convergence was examined. The 
observed data were fitted with the linear combination of the subd- 
warf model and a secondary template. At the beginning, a 50 per 
cent flux contribution was assumed for F and early G type compan- 
ions and the interpolation step size was 0.1 times the flux difference 
between spectra. At every decomposition iteration the interpolation 
step size was halved and after the fifth iteration it was set to 0.001. 
To avoid trapping, all template spectra were compared to the ob- 
served data at every iteration, but interpolation was performed only 
between the best-fitting template and its neighbours. 

The decomposition of late G and K companions required a 
more elaborated method. Their most significant effect on the flux 
distribution is a slight flattening of the flux in the red part of the 
spectrum and relatively strong Na D lines. In this case, the ini- 
tial flux ratio of the secondary was set to 15 per cent at the long 
wavelength end of the spectrum. Fitting started with a pre-defined 
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Figure 2. Panel a shows the V—J vs. J-H colour-colour diagram for our 
subdwarf sample. Optical magnitudes were collected from GSC2.3.2 and 
infrared from the 2MASS data base using VizieR. Apparently single stars 
aggregate near V— Jm— 0.5 and J—Hx—0. 1, while composite-spectra bina- 
ries show infrared excess (lower right corner). About 19 per cent of the stars 
have V—J > and J-H > 0. Open circles show synthetic colour indices of 
the 29 binaries resolved in this work. Typical errors of photometric colour 
indices are shown in the upper left corner. Panel b and c show the correla- 
tions of synthetic and observed V—J and J — H colour indices, respectively. 



section size to consider only spectral lines in the beginning, and 
the continuum was progressively considered by gradually increas- 
ing the section size up to three times the starting value. If the sec- 
ondary contribution decreased below 10 per cent, the spectrum was 
sampled with the smallest possible sections, which was set to three 
times the size of a resolution element. This strategy is necessary 
only for late-type companions and does not change the results for 
F and early G stars, which can be modelled based on spectral lines. 
The fitting windows were reset to their input configuration in the 
last step and kept fixed during error calculations, because chang- 
ing the size of fitting sections can introduce systematic shifts. In 
the case of late companions a high SNR ( > 100), medium to high 
resolution data, and good fluxing are critical. 

We consider that our semi-empirical approach in characteriz- 
ing the cool star is not just faster, but less ambiguous than work- 
ing with synthetic spectra for both components. To check if our 
interpolation returns real stellar spectra, the final templates were 
always compared to those calculated with the same parameters by 
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MILES Webtools 3 (Vazdekis et al. 2003). Companion spectral types 
were determined using the MILES classifications. We show an ex- 
ample for our fit and spectral decomposition of the sdB-F7V binary 
J2124+1506in Figure 3. 

To resolve double-lined close binaries of similar luminosities, 
individual radial velocities would be helpful. However, such a study 
would require high-resolution spectra and templates or calculated 
spectra for both components. A productive method to resolve these 
binaries would be to obtain both ultraviolet and optical spectra. 

3.4 Spectral fitting 

We used the same initial model with: T c ff = 40 000 K, logg = 5.6 
cm s~ 2 , [He/H]. = -1 and [(C,N,0)/H]. = -2 to fit our observa- 
tions. Maximum relative changes per iterations were limited to ±5 
per cent for T e g, ±2 per cent for logg and +100 per cent, -50 per 
cent for all abundances in order to maintain a stable TLUSTY con- 
vergence. We fitted only fluxed spectra with 80 A segments (each 
having six to 16 resolution elements depending on the resolution) 
using the entire spectral range. In order to decrease model atmo- 
sphere calculation time we used 30 depth points. Our tests showed 
this simplification affects model convergence before it would sig- 
nificantly change the emergent flux at low resolution. Model atmo- 
spheres were calculated in non-LTE radiative equilibrium without 
convection. Detailed profiles of H and He lines were included, but 
rotational broadening was not. A high rotational velocity is not typ- 
ical for subdwarf stars. Fitting of ~200 observed spectra took ~800 
h with six processors (15.2 GHz total) and needed the calculation 
of about 25 000 model atmospheres. 
All fits are available in our online catalogue 4 . 



3.4.1 Known issues 

The piecewise normalization procedure implemented in XTGRID 
can reduce the effects of continuum inconsistencies between the 
theoretical and observed spectra as well as the effect of interstel- 
lar reddening. However, the method creates discontinuities in the 
fit residuals of stars that suffer from large reddening. Such frac- 
tures in the residuals may divert the fitting procedure and intro- 
duce systematic shifts in the derived parameters. Therefore, we 
performed extensive tests to discover if a wavelength dependent 
extinction correction would change atmospheric parameters signif- 
icantly. We included the extinction function from Cardelli, Clayton 
& Mathis (1989) in XTGRID, and implemented an iterative redden- 
ing determination in the program. After recalculating some heav- 
ily reddened spectra, we did not find systematic differences. Our 
quoted error bars are larger than the effects of our normalisation on 
the parameters. We estimated the interstellar extinction from the fi- 
nal fits and listed the associated E(B-V) colour excesses in Table 2 
and 3. These values are consistent with the dust maps of Schlegel, 
Finkbeiner & Davis (1998) that provide an upper limit to the ex- 
tinction coefficient toward each target. 

Our binary decomposition method tends to attribute interstel- 
lar absorption lines along with other spectral lines to the template 
and fits the primary in this environment. This can cause degener- 
acy for binaries with similar components, like in 10716+2319. A 
single star fit provided a subluminous star with T e g = 12940 ± 300 

3 http: //www. iac . es/proyecto/miles/pages/webtools.php 

4 http: //pleione . asu. cas . cz/-nemeth/work/galex/ 
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Figure 3. Binary decomposition of J2124+1506. The single model (top panel) can not reproduce the composite spectrum, the G-band (~4300 A), Can H&K 
(3933 and 3968 A) and Mgi (5183 A) lines are clearly in the residual. The jagged model also indicates that the theoretical continuum is inconsistent with the 
observation. When the spectrum is modelled with an additional cool MS component (bottom), all major spectral features are reproduced well. The subdwarf 
temperature in this case was adjusted mildly, while there was significant change in the surface gravity and He abundance. A CNO abundance could not be 
measured in the spectrum. 



K and logg = 5.34*J - jji and slightly subsolar abundances. The ap- 
parent Ca and Mg lines, the H line profiles, the very asymmetric 
error in gravity and the infrared excess suggested the presence of a 
companion. Binary decomposition provided a much better fit with 
a B-Al V pair with r eff = 1 1 140 ± 180 K, log g = 4.391°* primary 
and T eff = 9310 K, logg = 3.67 secondary. 

Bright companions, like A stars, outshine sdB stars in the 
optical. Atmospheric parameters of relatively cooler subdwarfs 
in such binaries can be measured only with large uncertainty. In 
J2349+4119 the T eff = 7940 K, logg = 3.44, A4V companion 
contributes 87 per cent of light in the B photometric band. The de- 
termination of r cff = 11730 K and logg = 5.8 for the primary is 
less certain. Radial velocity measurements would help finding the 
true nature of this binary. 

In unresolved composite spectra XTGRID can attribute subd- 
warf metal lines to the lines of the companion, leading to incorrect 
CNO abundances. While C and/or N is relatively abundant com- 
pared to O and show lines in the spectrum, this affects the O abun- 
dance determination in the first place. Companion G-band absorp- 
tion overlaps with the strongest On lines (/14273-/14322 A) of cool 
subdwarfs in composite spectra binaries. These lines are important 
in the O abundance determination and their attribution to secondary 



features can, in some case, lead to an underestimated O abundance 
of subdwarfs. Contrarily, XTGRID might overestimate the O abun- 
dance of sdO stars in unresolved composite spectra binaries. Our 
current O abundances are only preliminary determinations, a com- 
plete study requires higher-dispersion spectra. 

The puzzle of the Nm/Civ /14640-/14665 A blend still chal- 
lenges the present spectral modelling. We observed the blend in 
about half of the hot He-sdO stars over a wide range of tempera- 
ture and our fits do not reproduce it satisfactorily. We revised our 
line list and updated the atomic data for the Civ /14658.3 A mul- 
tiplet with lines between the ls 2 5f(2F°)-ls 2 6g(2G), ls 2 5g(2G)- 
ls 2 6h(2H°) and ls 2 5g(2G)-ls 2 6f(2F°) levels from the Tubingen 
NLTE Model-Atmosphere Package (TMAP; Werner & Dreizler 
1999) available at the German Astrophysical Virtual Observatory 
page 5 . We noted a line strength inconsistency of the C m /14070 A 
multiplet with the C iv /14658.2 A multiplet as well as with the C iv 
/15801.3 and /15812 lines. Our models show too strong C m and too 
weak C iv lines compared to the observations, indicating an under- 
estimated temperature or an overestimated gravity. The effects of a 



5 http: //astro. uni- tuebingen. de/~TMAD/ 
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Figure 4. T e g~\ogg diagram. The ZAEHB and TAEHB (Dorman, Rood & O'Connell 1993) are marked with "a" and "b", respectively. The HeMS is taken 
from Paczyriski (1971) and Divine (1965) and are marked with labels "c" and "d", respectively. Stellar masses are marked along the HeMS in M Q . The grey 
lines are theoretical evolutionary tracks from Dorman, Rood & O'Connell (1993) for stellar masses from top to bottom: 0.480, 0.475, 0.473 and 0.471 M G . 
Line widths are proportional to evolutionary time-scales. Typical error bars are shown in the upper left corner. He-rich ([He/H], > -1) stars are indicated 
with filled circles, He-poor ([He/H], < -2.2) stars with filled triangles and stars with —2.2 < [He/H], < -1 are with filled squares. Subdwarfs in resolved 
composite spectra binaries are shown with open symbols following the same classification. Iso-luminosity curves for L/Ledj = -1.6, -2.0, -2.6, -3.0 and- 3.4 
are shown with dash-dotted lines. These are characteristic values in the luminosity distribution. 



change in gravity are low on metal lines, hence only a change in 
temperature could provide better fits. However, the ionisation bal- 
ance shows that the C iv ion is dominant over C in above ~42 000 K, 
hence an increase in the temperature would decrease the C in line 
strength, but would not increase the C iv opacity. Such changes in 
temperature or gravity would be inconsistent with the H and He line 
strengths as well. The high C and N abundance of these stars im- 
plies a high probability of heavier elements, in particular the iron- 
peak elements, which can be accounted for the missing opacities. 
Haas et al. (1996), Werner et al. (1998) and recently O'Toole & 
Heber (2006a) found high abundances of heavy elements in hot 
subdwarf stars. A review of the atomic data input is necessary and 
deserves a separate study. 



4 PROPERTIES OF HOT SUBDWARFS 

Our sample is large enough and appropriate for revisiting the dis- 
tribution of stars in the T cft - log g and T cft - He planes as well as 
to look for possible correlations in the He and CNO abundances. 
Such diagrams are very important tools in tracking subdwarf evo- 
lution and pulsational modes, because with the help of a model at- 
mosphere program, stellar temperature, gravity and element abun- 
dances can be derived without any assumptions on distances or stel- 
lar luminosities. A comparison of the distribution of stars in these 
diagrams is also useful to test our fitting method. 



Our spectral classification follows the simple scheme of Paper 
I: Subdwarfs with T sff < 40 000 K showing dominant H lines and 
weak He i are classified as sdBs; with T e ff > 40000 K and weak 
lines of He n are sdOs. Stars with dominant He i lines are He-sdBs, 
while those with dominant Hen lines are classified as He-sdOs. 
Stars with r eff < 20 000 K and log g < 5 are listed as BHB stars 
and with log g < 4.5 are classified as MS B stars. Using this cat- 
egorization, out of the 180 stars we found 166 subdwarfs, 124 of 
which are classified as sdB and 42 as sdO stars. The sdB/sdO num- 
ber ratio of ~3 is the same figure as previously determined (Heber 
2009). We found five He-sdB stars (or ~4 per cent) among the 124 
sdB and 26 He-sdO (or ~62 per cent) out of the 42 sdO stars. In the 
29 resolved spectroscopic binaries only four (namely: J07 10+2333, 
J1602+0725, J2020+0704 and J2038-2657) have sdO primaries 
and two binaries (J0047+0337 and J2331+2815) are in the sdB/sdO 
transition region. No He-sdO subdwarfs were found in composite 
spectra binaries. One system (J07 16+23 19) is a MS binary and one 
(J2349+4119) is a possible low mass pre-WD-MS binary. The re- 
maining 21 composites are classified as sdB-MS binaries. 



4.1 The T cff - log g plane 

The properties of our hot subdwarfs in the temperature-gravity 
plane are shown in Figure 4. The distribution of sdB stars follow 
the theoretical EHB and He-sdO stars are scattered near the He 
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main-sequence (HeMS), with some stars along post-EHB evolu- 
tionary tracks. The post-EHB tracks of Dorman, Rood & O'Connell 
(1993) for subdwarf masses from top to bottom: 0.480, 0.475, 0.473 
and 0.471 M G are shown in grey. The zero-age EHB (ZAEHB) 
and terminal-age EHB (TAEHB) are marked with "a" and "b", re- 
spectively. The HeMS is taken from Paczyhski (1971) and Divine 
(1965) and are labelled with "c" and "d", respectively. Stars with 
[He/H]. > -1 are marked with filled circles; most of these spread 
along the HeMS. Stars with [He/H], < -2.2 are indicated with 
filled triangles and gather at lower temperatures and gravities than 
stars with 2.2 < [He/H], < -1 that are shown with filled rectan- 
gles and crowd around T eS « 35 000 K and logg = 5.8. Resolved 
composite spectra binaries are marked with open symbols using 
the same scheme. Dash-dotted lines show iso-luminosity fractions, 
where the luminosity is expressed as the fraction of the Eddington- 
luminosity calculated from the effective temperature and surface 
gravity: 

log-^- =4xlogr cff - logg- 15.118. (1) 

J-Edd 

This distance independent quantity statistically separate stars. As- 
suming a constant stellar mass for sdB stars, these curves are also 
proportional to their absolute luminosities. 

Figure 4 reveals five distinct regions in the r c ff-logg diagram. 
The cooler, He-poor sdB stars (no. 1) are found around 28 000 K 
and logg = 5.45. These are potential slow pulsators {P « 40 - 170 
min). The hotter sdB stars (no. 2) near 33 500 K and logg = 5.8 are 
on average 10 times more He abundant and possible rapid pulsators 
(P as 2 - 6 min). The hot and He-rich sdO stars (no. 3) are found 
mostly between 40000 and 55 000 K near logg = 5.9 ± 0.4. Al- 
though numerous He-sdO stars are above the theoretical HeMS, as 
theoretical models predict, an asymmetric scatter toward high grav- 
ities is apparent in our data and noted in the literature (0stensen 
2009). He-weak sdO stars (no. 4) are located above both the EHB 
and HeMS. Finally, the least defined population: possible progen- 
itors of low-mass WDs, pre-WDs, post-Asymptotic Giant Branch 
(post-AGB) stars and central stars of planetary nebulae (CSPN), 
that might cross this region evolving towards their cooling tracks 
(no. 5) are also distinguished. The lack of stars along sections of 
fast evolution, in particular between -2.0 < L/L Edd < -1.6 is also 
notable in Figure 4. Furthermore, it is immediately seen in the fig- 
ure that He abundance, surface temperature and gravity are corre- 
lated. 

While sdB stars are assumed to evolve more or less along 
the indicated evolutional tracks, He-sdO stars in lack of a H enve- 
lope are expected to leave the HeMS towards higher temperatures 
and gravities after He exhaustion. However, this fast evolution can- 
not explain the observed clustering of high gravity He-sdO stars. 
Hu et al. (2008) found that stars with M ZAMS > 2 M Q can ignite 
core He-buming in non-degenerate conditions and after the EHB 
such stars can develop a He-burning shell and a different struc- 
ture than canonical subdwarfs. An observational property of these 
stars would be their clustering below the theoretical HeMS. From 
an observational point of view, incomplete line broadening param- 
eters of He lines can also be responsible for overestimating gravity 
in such stars. More observations are needed to better characterize 
these high gravity He-sdO stars. 

A known shift in sdB temperatures by about +2000 K com- 
pared to the blue edge of the theoretical EHB band (Charpinet et 
al. 2007) is also noticeable in Figure 4. Similar shifts have been ob- 
served in previous analyses with various modelling and fitting tech- 
niques in volume limited samples. Hence, this discrepancy deserves 
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Figure 5. Subdwarfs in the theoretical instability region. Contours were ren- 
dered from Charpinet et al. (2011) and are overplotted on our data. Rapid 
(p-mode) pulsators are concentrated within the subdwarf instability region 
(Charpinet, Fontaine & Brassard 2001), indicated by the contours, driving 
is most efficient in the three highest contours (solid lines). The dash-dotted 
line separates slow and rapid pulsators. This boundary was found observa- 
tionally based on an independent dataset (Charpinet et al. 2010). 



attention from the theoretical side. Jeffery & Saio (2006) showed 
that including more iron-peak elements in the models would shift 
the theoretical EHB towards the observations. 

The contours in Figure 5 mark the theoretical pulsational in- 
stability region by Charpinet, Fontaine & Brassard (2001) and are 
overplotted on our sample. The pulsations are driven by the iron- 
peak opacity bump. Charpinet, Fontaine & Brassard (2009) and Hu 
H. (2011) calculated models with equilibrium Fe abundance be- 
tween radiative levitation and gravitational settling. The contours in 
Figure 5 were reproduced from Charpinet et al. (201 1) and show the 
number of pulsating £ = 0, p-modes. The outermost contour corre- 
sponds to a single mode and the highest in the centre to seven. Driv- 
ing is most efficient in the three innermost contours. Possible rapid 
(p-mode) pulsators are located around 33 500 K and logg = 5.7 
while presumed slow (g-mode) pulsators are around 28 000 K and 
logg = 5.4. The boundary that separates slow and rapid pulsators 
is marked with the dash-dotted line and was found independently 
(Charpinet et al. 2010). We refer to these groups as possible pul- 
sators, although variable and non- variable stars co-exist in their ob- 
served parameter space. 0stensen et al. (2010) found a ~10 per 
cent incidence rate for rapid pulsators in their survey for pulsating 
hot subdwarfs with the Nordic Optical Telescope. Slow pulsators 
show a larger incidence rate, reaching ~75 per cent according to 
Green et al. (2003), which was recently confirmed by the Kepler 
sdB sample (Ostensen et al. 2011). In between these two groups 
are the hybrid (or mixed-mode) pulsators that show both slow and 
rapid modes. Based on the distribution of stars in Figure 5 many of 
our subdwarfs are possible pulsating candidates and a deep astero- 
seismic follow-up would be reasonable. Based on the locations of 
known pulsators we sorted out and flagged candidates for pulsation 
studies in Table 2. We would like to note that the GALEX sample 
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Figure 6. Left: He abundance versus effective temperature. Two trends of increasing abundance with effective temperature can be seen. Also remarkable are 
the clumping of possible slow (#1, near 28 000 K and [He/H], = -1.5) and rapid sdB pulsators (#2, near 33 500 K and [He/H], = -2.7) and the hot He-rich 
sdOs (#3). The sequence of He-weak subdwarfs (#4 and dash-dotted line) follow a similar correlation with effective temperature as more He rich stars (short 
dashed). The short full line is the best fit from Edelmann et al. (2003) for the He-weak sequence. Suspected progenitors of low-mass WDs and BHB stars are 
also possible in this region of the HRD (#5). The pure He atmosphere of J085 1-1712 is an extreme case, we determined only a lower limit of its He abundance. 
Upper limits are indicated with open down triangles and binary stars are encircled. Right: Same data without error bars and point sizes are proportional to 
surface gravities (point size = log g - 4.8) as shown in the upper left corner. Indicated regions are intended to guide the eye and are not derived from the data. 
Interesting to note the distribution of composite spectra binaries as well as surface gravity correlations. 



presented here is a consistent set, the relative position of the stars 
to each other are free of major systematic shifts. 



4.2 The T cff - He diagram 

The He abundance in hot subdwarf stars is another fundamental 
parameter and must be investigated along with surface temperature 
and gravity. The panels of Figure 6 show the distribution of stars 
in the T cff - He diagram and their correlations with surface gravi- 
ties. The five groups introduced in Section 4. 1 are separated in the 
Teg-He plane as well. The He-poor (no. 1) and He-rich (no. 2) sdB 
stars, the sequences of He-sdO (no. 3) and He-weak sdO (no. 4) 
stars are marked in Figure 6 as well as BHB stars and possible WD 
progenitors (no. 5). The two r c ff-He correlation sequences first de- 
scribed by Edelmann et al. (2003) are also remarkable. In general, 
the He abundance increases with temperature. Our He detection 
limit is at [He/H], as -3 and does not allow to trace the two helium 
sequences deeper, although some data suggest a continuation of 
the He-weak sequence. Similarly large and homogeneously mod- 
elled samples like the Hamburg Quasar Survey (HS, Edelmann et 
al. 2003), Supernova la Progenitor surveY (SPY, Lisker et al. 2005, 
Stroer et al. 2005) and Sloan Digital Sky Survey (SDSS, Hirsch, 
Heber & O'Toole 2008) confirm this observation (O'Toole 2008). 
In Figure 6 we plot the best fit line (short dashed) for sdB stars from 
the HS sample (Edelmann et al. 2003) for the He-rich sequence: 

[He/H]. = -3.53 + 1.35 (yj^-2.0()J. (2) 

This independent trend fits our sdB data as well, in particular for 
He-rich sdBs and passes through the region of He-sdO stars. We 
also plot their regression (full line) for the He-weak sequence. 



However, we found a different trend (dash-dotted line) for these 
stars, possibly because the temperature range of the He-weak se- 
quence was under-represented in the HS survey (T e g < 40000 K). 
Separating He-weak stars at 

[He/H]. < 1.08 x 10- 4 r eff - 6.64, (3) 

and fitting them in a similar manner we got: 

[He/H]. = -4.26 + 0.69 1^ - 2.0oJ . (4) 

Although, it can be seen from these fits that the distribution of sub- 
dwarf stars is much more complex and cannot be satisfactorily de- 
scribed with simple linear trends, hence we omit their significance 
test. Multiple breaks can be found in the He abundance distribu- 
tion. There is a clear cut in temperature (vertical dashed line) in 
Figure 6, separating sdB and sdO stars somewhere between 36 000 
and 40 000 K depending on the He abundance. To the right of this 
line are the two, relatively well defined groups of sdB stars, while 
to the left are the hot, either He-weak or He-rich sdO stars. A He 
abundance split is also remarkable at the line [He/H], = -0.5 + 0.5 
(horizontal dashed line), in particular over 40 000 K, where no stars 
can be found. Such a gap can be seen in the SPY data as well. Ob- 
servationally this corresponds to a region where stars are expected 
to either never enter or pass through quickly during their evolution. 
The major formation theories in Section 7 also predict a lack of 
stars in this region. Similarly to the T eff - log g distribution, possible 
rapid and slow sdB variable stars are separated according to their 
He abundance. Rapid pulsator candidates with [He/H], as —1.5 
have, on average, an order of magnitude higher He abundance then 
their slow pulsating counterparts with [He/H], as -2.7. A region 
of photospheric mixing and He enrichment near T e g as 38 000 K at 
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Figure 7. Sample fits (grey lines) of five subdwarf spectra (black lines) with increasing effective temperatures from bottom to top. The strongest C and N lines 
are marked. 



logg « 6 was predicted by Groth, Kudritzki & Heber (1985) and 
recently reviewed by Miller Bertolami et al. (2008). Our data show 
some abundance extremes in this region with [He/H], > 1 for He- 
sdO stars, culminating at the location of J085 1- 17 12 in which we 
could determine only a lower limit of the He abundance. Toward 
higher temperatures from the sdB/sdO border the He abundance 
in He-sdO stars shows a decrease with temperature. The He abun- 
dance in these stars shows a correlation with mass along the HeMS. 

Below the sdB/sdO border, a few He-sdB stars in our data 
show a He overabundance, but in general we can conclude that the 
He abundance is proportional to the effective temperature in sdB 
stars, and no stars can be found with similarly high He abundances 
like in He-sdO stars. In the proximity of the sdB/sdO transition both 
He-rich and He-poor sdB and sdO stars can be found. The case of 
He-weak sdO stars is simpler. They show a correlation with effec- 
tive temperature similar to sdB stars and are connected to sdB stars 
at the sdB/sdO border. Their binary frequency and surface gravity 
diversity suggests that the He-weak sequence is a mixture of mul- 
tiple populations. We discuss this point further in Section 7. 

While parameter errors of temperature and gravity show grad- 
ual increase to higher temperatures, the error budget for the He 
abundance seems to be three-fold: sdB and BHB stars have well de- 
fined atmospheric parameters with symmetric error bars. He- weak 
sdO stars have significantly larger errors and show asymmetries in 
temperature. He-sdO stars show a rather symmetric temperature, 
but very asymmetric He abundance errors; our method tends to 
stop at the lower limit. Probably because at such high temperatures 
and He abundances the spectrum is He dominated and hydrogen 
lines are difficult to trace with low-resolution spectroscopy. In such 
conditions He becomes the reference element and the x 1 is l ess 



sensitive for the change of the shallow He lines. Beyond the low 
resolution of our spectra, the incomplete chemical composition of 
our models can also be blamed for the larger parameter errors of 
He-sdO stars. 

By comparing the panels in Figure 6 further differences be- 
tween the two He sequences become apparent: there is a strong cor- 
relation between surface temperature, gravity and He abundance of 
sdB and He-sdO stars, but this cannot be said about He-weak sdO 
stars. Among these stars a much weaker correlation and scattered 
high gravity stars can be found. Also notable are the low tempera- 
ture sdB and BHB stars that have a wide distribution in He abun- 
dance and low surface gravities. Some of our subdwarfs show a 
low He abundance at a relatively high gravity, these are possible 
progenitors of low-mass WDs. 



4.3 CNO trends 

Metal abundance patterns are difficult to obtain from optical spec- 
tra, because the vast majority of CNO lines are in the ultraviolet 
and CNO lines in the optical are relatively insignificant. In spite of 
the weakness of these lines, our analysis provided C and N abun- 
dances for a number of sdB and sdO stars. Our metal abundance 
determination is based on the global effects of a specific chemi- 
cal composition on the entire spectrum as well as on individual 
lines. We adopted this method mainly because in low-dispersion 
optical spectroscopy weak metal lines are heavily blended and hard 
to identify, but, thanks to the flux blocking of UV lines the over- 
all shape of the energy distribution is also affected by the chemical 
composition. Sample fits for various subdwarf types are shown in 
Figure 7. Following the C&N analysis of sdO stars by Stroeer et al. 
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(2007), we summarise our results for the GALEX sample in Fig- 
ure 8. In the r eff - logg diagram the coolest EHB and BHB stars 
(T e fF < 24000 K) show both C and N abundances. Cooler sdB stars 
between 24 000 K and 30 000 K show a higher N abundance, but 
no C in general. In a narrow range around 30 000 K both C and N 
were found in some stars. The majority of sdB stars over 30000 K 
show much less diversity, in most of these stars only He was found 
similarly to post-EHB stars. He-sdO stars show a diversity of sur- 
face abundances, but a pattern can be outlined. N-class stars are 
concentrated below r e g < 45 000 K. He-sdO stars with higher tem- 
perature and surface gravity belong to the C-class, while a group of 
C&N-class objects can be found in a narrow strip at ~45 000 K. The 
C-class seems to belong exclusively to hotter He-sdO stars. Finally, 
the hottest sdO stars (T s g > 60000 K) show a N overabundance in 
their spectra. The evolutionary status of these stars is an interesting 
question. They might belong to evolved He-sdO stars moving away 
from the HeMS as well as to He-weak sdO stars moving towards 
the HeMS in the temperature-gravity diagram (Husfeld et al. 1989). 

Based on the T c g - He graph, we can confirm the conclusions 
of Stroeer et al. (2007) that He-weak sdO stars do not show metal 
traces. We found only a few stars around T eff as 40000 K near 
[He/H], a -0.5 where a small clumping was observed in the SPY 
sample. However, in our data, the trend of decreasing He abun- 
dance with increasing effective temperature can be traced down to 
[He/H], > in contrast with [He/H], > 1 found by Stroeer et al. 
(2007). 

Figure 8 gives only a qualitative picture of the C and N distri- 
bution, because these graphs do not reflect element abundances. To 
further investigate the C and N abundance distributions and possi- 
ble C-He, N-He and C-N abundance correlations we refer to the 
panels of Figure 9, where only positive determinations are shown. 
In Appendix A we repeat these figures supplemented with error 
bars and upper limits. 

Panels a, b and c of Figure 9 show CNO trends with respect to 
temperature. The approximate border between sdB and sdO stars at 
T c ff ~ 38 000 is shown by a vertical line and solar abundances from 
Grevesse & Sauval (1998) are indicated with horizontal dashed 
lines. Out of the 180 stars, we measured a He abundance in 149; 
68 showed N, 39 showed C, but only 7 showed O in their spectra. 
A quick look reveals that the bulk of sdB stars have subsolar C and 
O abundances and roughly solar N abundances. He-sdO stars show 
higher abundances, strongly correlated with effective temperature. 
Interpreting the abundance distribution in the direction of increas- 
ing temperature: both C and N show a similar abundance in cooler 
sdB or BHB stars in the range 20 000 K < T slf < 23 000 K of about 
[C/H]., [N/H]. ~ -3. In sdB stars between 23 000 K < T cff < 
34000 K a lower abundance of about [C/H]., [N/H]. « -4 was 
measured. Above r cff « 35 000 K across the sdB/sdO border a sig- 
nificant abundance increase of C and N is observed reaching a peak 
near r eff a 43 000 K at [C/H]., [N/H]. a -1, over this tempera- 
ture the abundance of both elements decrease steadily with temper- 
ature. The hottest stars indicate a gradual abundance increase over 
r c fj a 60 000 K. Only a few stars without He in their spectra show 
C or N. These trends are very similar to those found in the He abun- 
dance distribution and similar observations were made in the SPY 
data (Lisker et al. 2005, Stroeer et al. 2007). 

In panels d and e of Figure 9 the C and N abundances are 
plotted with respect to the He abundance. The 39 stars that have 
C determination have He as well. Out of the 68 stars showing N, 
64 show a He abundance. The stars in which both C and N were 
found simultaneously are fitted with a full line. A similarly com- 
plex C and N abundance pattern can be outlined from this data as 
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Figure 8. Metal abundances in the T c ff - log g and T^g - He diagrams. 



in the case of He. In sdB stars both C and N show large variations, 
although a higher N abundance in cooler sdB stars is apparent. To- 
wards higher He abundances, starting at [He/H]. a — 1, He-sdB 
stars nicely connect sdB and He-sdO stars (Ahmad & Jeffery 2003) 
and show a higher C and N abundance than sdB stars. Unlike in sdB 
stars, the He abundance in He-sdO stars decreases with increasing 
effective temperature. Therefore, in the direction of increasing He 
abundances, we observe the hottest He-sdO stars first, and finally, 
we see the most He abundant and cooler He-sdO stars. We observed 
a similar trend for the C and N abundances as well. Due to this dif- 
ferent population, an abundance drop is observed at [He/H], a 0, 
which is followed by a steady abundance increase with increasing 
He abundance. These trends reach a peak at [C/H]., [N/H], a -1 
and were fitted with a dashed line in panel d for the C abundance: 

[C/He]. = 1.72(±0.23) [He/H]. - 3.55(+0.17) (5) 

and with a dash-dotted line in panel e for the N abundance: 

[N/He]. = 1.00(±0.13) [He/H]. - 3.16(±0.18), (6) 

respectively. These fits suggest that He-sdO stars with a low He 
abundance are either C-rich or N-rich. Over [He/H], a 0.55 both 
C and N is present in the spectra and C is more abundant. Finally, 
the most He abundant He-sdO stars show only N in their spectra. 
The hottest sdO stars complicate the trend, but it can be seen that to- 
wards higher He abundances (lower temperature) the N abundance 
progressively remains below the C abundance. The ratio of these 
two elements has important connections with formation theories. 
In the case of the extremely He-rich star J0851-1712 we could not 
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Figure 9. Panels a, b and c show the measured C, N and O abundances with respect to effective temperature. A vertical line at 38 000 K separating sdB and 
sdO stars is overplotted as well as the solar abundance for each element with a horizontal line. Panels d and e show correlations between the C, N and He 
abundances. In panel /, the C-N abundance correlations are shown, trends are fitted for N-rich stars (full line), C-rich stars (dashed line) and both C- and N-rich 
stars (dash-dotted line). 
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determine an upper limit to the He abundance to support this trend. 
For such He dominated stars and He-sdO stars in general, it might 
worth considering changing the reference element to He because 
abundances relative to H may show extreme values merely due to 
the vanishing H in the atmosphere. 

In panel / the N abundance is plotted with respect to the C 
abundance. The three lines show the linear fits to stars that show a 
C (dashed line), N (full line) and both C and N abundances (dash- 
dotted line). While data for C or N-rich stars independently show 
about the same abundances in He-sdO stars, stars near the C-N 
peak that show both C and N simultaneously suggest a C to N ratio 
of about ~3. Considering the asymmetric errors of our abundance 
determination (see Appendix A) this number is a lower limit. The 
same trend predicts a N overabundance in the coolest sdB stars. 
We found equilibrium abundance at [C/H]., [N/H]. a -3.5 that 
corresponds to [He/H], a -1.25; just above the group of possible 
short period sdB pulsators in the temperature-He-abundance dia- 
gram. C&N-class objects are C-rich above this He abundance and 
N-rich below. Based on the number of detections the C and N abun- 
dance peaks of He-sdO stars are also recognisable in panel /near 
[C/H], a -2.5 and [N/H]. a -2.3, respectively. We note that our 
C abundance determination in some cases has asymmetric errors 
and the N abundance has large errors making the derivation of fine 
details difficult. However, our observations confirm the C&N clas- 
sification scheme of He-sdO stars. Interestingly, the most He abun- 
dant star ( J085 1-1712) showed a N abundance, but we could not de- 
termine a C or O abundance. It might belong to the N-class. Zhang 
& Jeffery (2012) showed that WD mergers can produce such abun- 
dance patterns. Their slow mergers are predicted to be N-rich and 
populate the low-mass end of the HeMS around T^g a 41 000 K. 
In turn, fast mergers are expected to be C-rich with higher masses 
around r eff a 45 000 K. Their composite model retains both ele- 
ments with N overabundance in low-mass and C overabundance in 
high-mass He-sdO stars. Such C enrichment can be explained by 
the standard hot-flasher scenario as well, for He-sdO stars it pre- 
dicts a C/N ratio in the range of -0.85 to ~9.5 (Miller Bertolami et 
al. 2008) in agreement with our observations. 

In Appendix A we show upper limit determinations and er- 
ror bars. These graphs show even more similarities with the He 
abundance distribution, but a deeper look reveals some differences 
as well. The appearance of the C-weak and N-weak sequences is 
remarkable. The C-N peak shows a structure: C shows a peak at 
r efr a 46000 K and overabundant compared to N up to r eff a 
55 000 K, while N shows a peak at T^g a 42 000 K and decreases 
with increasing temperature faster than C. This structure has al- 
ready been observed in He-sdO stars (Hirsch 2009) and was in- 
vestigated in the case of double He- WD mergers (Zhang & Jeffery 
2012). From the model analysis of Hirsch (2009) we found a mean 
temperature of -44 600 K for C-rich and -40600 K for N-rich He- 
sdO stars, in agreement with our results. In the hottest He-sdO stars 
(7eff > 60 000 K) we were able to determine N and found a slight 
anti-correlation with the He abundance unlike at lower tempera- 
tures. Our O determination is not yet reliable to derive such trends, 
but a correlation with the He abundance is possible. 

The boundary that separates rapid and slow sdB pulsators 
(dash-dotted line in the top panel of Figure 8) correlates with 
the observed abundance patterns. Suspected rapid pulsators near 
33 500 K are He-rich and do not show metal traces, while slow pul- 
sators around 28 000 K are He-poor and many show a N abundance. 
We found 27 stars in the temperature-gravity region of known rapid 
pulsators. Out of the 27 stars, we found a N abundance in six (-22 
per cent). Similarly, we found 23 stars at the location of known 
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Figure 10. Same r c ff-logg diagram as Figure 4, here showing only com- 
posite spectra binaries. A crowding of composites can be observed near the 
location of possible rapid pulsators at T^g = 33 500 K, log g = 5.8. 

slow pulsators and 16 (-70 per cent) of these show a N abundance. 
In between these groups, where hybrid pulsators are expected, we 
found both C and N in most of the stars. Interestingly, flashless 
evolutional tracks by Hu et al. (2008) connect the group of C&N 
class sdB stars with C&N class high gravity He-sdO stars. Recently, 
Geier et al. (2012) found a He isotope anomaly ( 3 He is strongly en- 
riched in the atmosphere) in eight stars in a sample of 46 sdB stars. 
These stars are situated in a narrow strip in the temperature-gravity 
plane between 27 000 and 31 000 K, where we see C&N class sdB 
stars in our survey. No connection has been found in earlier stud- 
ies between surface abundance patterns and pulsational instability 
(O'Toole & Heber 2006b). However, we observed different He, C, 
and N contents along the EHB that might play a role and need fur- 
ther investigations. 

4.4 Binaries 

Figure 10 shows the distribution of subdwarfs in composite spec- 
tra binaries in the temperature-gravity plane. We did not find He- 
sdO stars in such binaries. Based on the Ca infrared triplet Aznar 
Cuadrado & Jeffery (2001) and Aznar Cuadrado & Jeffery (2002) 
found that companions from composite spectra are in the range of 
0.8-1.2 M Q F-G type MS or subgiant stars. We confirm this obser- 
vation and list the results of our binary decomposition in Table 4. 
Monochromatic flux ratios are listed in Table 5 for the companion 
stars at the effective wavelengths of the broadband Bessell filters 
(Bessell 2005). In parentheses we also give predicted flux ratios 
that were derived from the synthetic subdwarf spectrum and com- 
panion template outside our observed spectral range. 

Before starting spectral decomposition we always performed 
a single star fit and used the results in the input of the binary model. 
This allowed us to estimate the shifts in subdwarf atmospheric pa- 
rameters during decomposition with respect to a single star fit. We 
emphasise that the relative shifts in Table 4 are specific to our mod- 
elling and fitting method that does not make assumptions for the 
companion types. Temperature shifts show large variations with an 
average of - 1300 K increase. Surface gravity shows a stronger cor- 
relation, with an average 0.7 dex increase. Similar strong correla- 
tion can be observed in the He abundance shifts. On average a 0.42 
dex increase with a small scatter is indicated by the spectral de- 
composition of the 14 stars in which He was measurable in both 
the single star and binary fits. Spectral decomposition helped to 
achieve a He abundance measurement in four cases while in six 
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Table 4. Parameter shifts during decomposition. He abundance determination from upper limit to positive identification (normal font), upper limit before-after 
(italic font) and positive identification before-after (bold font) are distinguished. Companion absolute V magnitudes are derived from spectral types and from 
binary decomposition for subdwarfs. Absolute magnitudes of binary components allowed to estimate distances. 
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cases it confirmed an upper limit measurement. We adopted abso- 
lute V magnitudes based on the spectral type and luminosity class 
of the companions from Gray (1992). Using the flux ratios from Ta- 
ble 5 we calculated subdwarf absolute V magnitudes and approxi- 
mate distances. The absolute magnitude distributions of subdwarfs 
and their companions are shown in Figure 1 1 and confirm a consis- 
tent spectral decomposition. F type MS stars are the easiest com- 
panions to find, hence these are the most frequent in our sample. 

We did not find line asymmetries indicating a measurable ra- 
dial velocity difference between the components in any of our stars. 
Radial velocity did not show variations in the much larger, low- 
resolution sample of similar composite spectra binaries of Green et 
al. (2008), which suggests that orbital periods must range from at 
least a few months to few years. Recently, 0stensen & Van Winckel 
(2012b) measured radial velocity variations and detected orbits of 
eight sdB-F/G binaries. They found orbital periods in the range of 
500-1200 days. 

The effects of spectral decomposition of subdwarfs with F and 
early G type companions were found to be a positive shift in the 
subdwarf temperature, gravity and He abundance. In the case of late 
G and early K-type companions, due to the lower flux contribution, 
we expect lower shifts. A shift to higher gravities, hence to lower 
luminosities of some of these sdB stars is possible. A flux flattening 
and strong Na D lines can also indicate binarity. The high SNR 
allowed us to detect signs of such companions in the far red part 



of the spectra. However, such faint stars do not contribute enough 
flux next to a 40-300 times more luminous subdwarf to carry out 
a proper decomposition. Due to the lack of strong H lines in their 
spectra and their low contribution, these stars are not expected to 
significantly affect the derived subdwarf atmospheric parameters. 
The binary population with late MS (K and M dwarf) companions 
is high according to radial velocity surveys and probably similarly 
high with WDs. Such companions are too faint to show resolvable 
spectral signatures. 

Because close binaries are usually found by radial velocity 
surveys we checked the absolute shifts of our spectra relative to 
synthetic spectra. Figure 12 shows the radial velocity distribution 
in the GALEX sample with respect to the kinematic local standard 
of rest (LSR) frame. We found a two component Gaussian velocity 
distribution with velocity dispersions cr { = 82 and cr 2 = 26 km s" 1 , 
and mean velocities vi = -25 and V2 = -3 km s" 1 , respectively. 
Half of the stars in the catalogue show a radial velocity vv(LSR) > 
47 km s~' and 42 show v r (LSR) > 100 km s" 1 which, in some case, 
may be evidence of binarity. These stars are distinguished with the 
"rv" flag in Table 2. Two subdwarfs (J0321 +4727 and J2349+3844, 
Kawka et al. 2010) were found in close binaries and J1717+6757 
is a double-degenerate binary (Vennes et al. 2011b) out of the 42 
kinematically peculiar stars. J0716+2319 shows binarity and will 
be discussed in a forthcoming paper. Interestingly, four out of the 
29 composites show a notable radial velocity. In fact, some stars can 
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Table 5. Monochromatic flux ratios (F^ c , 



J (Fa, 



+ F,i,id)) of binary 



companions at the effective wavelengths of the broadband Bessell filters. 
Flux ratios outside our observed spectral ranges were calculated from syn- 
thetic and template spectra, these predicted flux ratios are listed in paren- 
theses. 
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be hierarchical triples (Maxted et al. 2001), these require further 
observations. 

We found eight composite spectra binaries among the 27 stars 
in the empirical region of rapid pulsators, while only one binary out 
of the 23 stars in the region of slow pulsators. The higher frequency 
of composite spectra binaries among possible rapid pulsators is ob- 
vious and supports the importance of binary evolution in subdwarf 
formation. A possible correlation between binarity and rapid pul- 
sations requires further investigations. 

Figure 13 shows the SEDs of the components in the sdO- 
G3.5III binary J2038-2657. The sdO dominates the flux in the UV, 
the subgiant emerges in the optical and dominates the infrared. 
The inset shows our observation, the components were smoothed 
to the correct spectral resolution. Our recent, high-resolution time- 
resolved spectroscopic follow-up revealed variable H tt emission in 
this binary. 



5 LUMINOSITY AND MASS DISTRIBUTION 

The GALEX sample presented here is one of the largest homo- 
geneously analysed subdwarf sample to date, therefore we de- 
scribe the cumulative luminosity distribution function (CLDF) of 
the 180 stars shown by the thick histogram in Figure 14. The bot- 
tom axis is the Eddington luminosity-fraction calculated from the 
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Figure 11. Absolute V magnitude distributions of hot subdwarfs and their 
companions. 
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Figure 12. Velocity distribution relative to the kinematic LSR. 



derived effective temperatures and surface gravities. The shaded 
area shows the probability density function calculated from the 
smoothed CLDF. The luminosity is also indicated in solar units on 
the top axis, assuming an average stellar mass of 0.48 M Q . This 
scale is valid only for sdB stars and their descendant sdO stars, and 
represents a lower limit for He-sdO stars that have progressively 
higher masses toward higher temperatures along the HeMS (Di- 
vine 1965, Paczynski 1971). Again, we found five well separated 
sections in the luminosity distribution. We identify the lowest lu- 
minosity stars below log(L/LEdd) < -3.5 with possible progenitors 
of extremely low mass WDs such as J0805-1058. The group be- 
tween -3.5 < log(L/L E(W ) < -3 consists of possible BHB stars, 
like J0746+0610, J1738+2634 and J2153-7004. About 55 per cent 
of our sample have luminosities between -3 < log(L/LEdd) = -2.6, 
we identify these with most of the sdB and some of the low-mass 
He-sdO stars. At higher luminosities the slope of the CLDF is de- 
fined by sdO and He-sdO stars up to log(L/Ledd) = -2. Finally, 
the most luminous sdO, post-AGB, or CSPN stars are found at 
log(L/L Ed d) > -1.6. After passing their maximum light, these stars 
become WDs and cross the subdwarf region on their way to the 
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Figure 13. Spectral energy distributions of the components in J2038-2657. The TLUSTY SED for the sdO star (grey) and the G5III template from the HILIB 
library are marked. Bessell B, V and 2MASS J and H filter transmission curves cut out the contribution of each star to the synthetic photometry. Inset: 
Magnified part of the spectrum where decomposition was performed. The HILIB G5III template is in grey, overplotted with our interpolated MILES template 
in black, their match is excellent. The G5III and the sdO spectra adds up nicely and this binary model fits the observation. The binary model is shifted up for 
clarity. At the effective wavelengths of Bessell B, V and 2MASS J and H filters we show GSC2.3.2 (quick- V) and 2MASS J and H photometry with black 
diamonds. Data for the Bessell B filter is the associated optical source magnitude form 2MASS. Grey circles show our synthetic photometry. 
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Figure 14. Luminosity distribution of 1 80, mostly hot subdwarf stars from 
the GALEX survey. The locations of the ZAEHB, TAEHB and terminal 
age post-EHB (TAPEHB) are marked. Stellar luminosities are expressed in 
solar luminosity on the top axis assuming a subdwarf (sdB, sdO) mass of 
0.48 M Q . He-sdO stars show a wide mass distribution. 



lowest luminosities. Two peaks at higher luminosities are notice- 
able. We identify the peak at log(L/L E[ i d ) w -2.5 with He-weak 
post-EHB sdO stars and at log(L/LEdd) ~ -2.15 with He-sdO stars. 
These groups can be seen to separate in Figure 15, where anal- 
ogously to the HRD, the distribution of subdwarfs, WDs (Eisen- 
stein et al. 2006) and stars from the MILES library are shown in the 
Eddington-luminosity fraction versus temperature plane. The He 
content and binary frequency on the He-weak sdO sequence seems 



to be different for the two sdB groups and suggests that He-rich 
sdB stars reach higher temperatures in their post-EHB evolution 
than He-poor sdB stars. 

To estimate the mass-distribution in our sample we compared 
our observed CLDF to a synthetic CLDF calculated from the sub- 
dwarf evolutionary tracks of Dorman, Rood & O'Connell (1993). 
This method looks for the most probable mass-distribution to re- 
produce the CLDF similarly to the method of Zhang, Chen & Han 
(2010). The Dorman evolutionary tracks were calculated for eight 
core-masses at eight different surface He abundances with various 
envelope thicknesses, giving 137 tracks altogether. However, they 
do not support the entire observed temperature and gravity range. 
The tracks were calculated for the HB with the EHB at its hot 
extreme that, along with the HeMS, is only partially represented. 
For this reason, we fitted only stars with -3.2 < log(L/L Edd ) < 
- 1 .7 and excluded low temperature and high gravity He-sdO stars 
(Tuff < 45000 K and logg > 5.9). This way we also excluded 
possible low-mass pre-WDs in the sample that are not supported 
by the tracks. The CLDF of the remaining stars is shown by the 
thin histogram in Figure 14. Next, we trimmed the evolutional 
tracks to the observed temperature and gravity range. This step 
was necessary, because due to the approximately constant lumi- 
nosity along the HB and EHB, a wide range of temperatures and 
gravities can contribute to a given point in the CLDF, causing a 
degeneracy in our solutions. Then, from the effective temperature, 
surface gravity and the age from the zero-age horizontal branch 
(ZAHB) of the Dorman-models we calculated evolutionary tracks 
in the Eddington-luminosity fraction versus age plane. Next, we 
interpolated these tracks to a new and common age axis. Finally, 
we applied &x 2 minimisation to determine the mass-distribution in 
the sample using the normalised weighted average of these tracks 
(synthetic CLDF) and fitted the observed CLDF. Such an analysis 



also allows for estimating the evolutionary times on the EHB and 
post-EHB phases. 

We assumed that our sample is complete in the observed lumi- 
nosity range and well-distributed, representing the true frequencies 
of all subdwarf types. Our main goals were to describe the slopes 
and the break points of the CLDF in general. The models repro- 
duce both sdB and sdO stars reasonably well, although they pre- 
dict a slower evolution and higher luminosities for post-EHB stars. 
From the theoretical models we found an average EHB lifetime of 
~100 Myr, followed by ~50 Myr post-EHB evolution, in contrast 
with Zhang, Chen & Han (2010) who found 160 and 20 Myr, re- 
spectively. The observed CLDF also suggests that the post-EHB 
evolution is about 50 per cent of the EHB lifetime. However, large 
discrepancies are expected if the observed sample is not complete 
or not fully described by these models. More recently, 0stensen 
et al. (2012a) found -80 Myr EHB lifetime followed by -62 Myr 
post-EHB evolution in the case of the He-rich sdB pulsator KIC 
1718290. 

We found that our subdwarfs have a well defined mass distri- 
bution around M = 0.52 ± 0.02 M G . A small contribution of stars 
with M = 0.48 ± 0.01 M Q and a shallow distribution up to 0.6 M Q 
was also found. We would like to note that the -2000 K system- 
atic shift found for sdB stars comparing to these theoretical models 
would shift the CLDF for sdB stars by about -0.1 unit to lower 
luminosities. This adjustment would shift the mass distribution of 
sdB stars to considerably lower masses by about 0.1 M . 

By applying theoretical models to the SPY sdB and sdO 
data Zhang, Chen & Han (2010) derived the number of sdO stars 
evolved form sdB stars to be observationally 12 per cent. Using our 
spectral classification criteria and the class separation in the tem- 
perature He abundance diagram (Figure 6) we found 120 sdB and 
15 sdO stars, or about 12.5 per cent that are possible post-EHB 
stars. 



6 MODELLING HOMOGENEITY 

Systematic shifts inevitably arise when data from different instru- 
ments are modelled with various model atmosphere codes. To de- 
crease such shifts in the atmospheric parameters of subdwarfs some 
aspects need consideration, such as: various modelling assumptions 
and implementations, different atomic data, various fitting meth- 
ods and error estimations can all introduce systematic differences. 
Therefore, independent model atmosphere analyses on large and 
homogeneous datasets, like the HS, SPY and SDSS data would be 
desirable with various modelling and fitting codes to find system- 
atic offsets and their reasons. 

Beyond homogeneous modelling and fitting techniques, con- 
sistent atmospheric parameters would require homogeneous data as 
well (similar SNR, resolution, spectral range). This requirement is 
difficult to comply due to observing time restrictions. For this rea- 
son a careful reduction is important, in particular, to apply a correct 
dispersion solution because even small departures from linearity or 
an incorrect radial velocity correction can change the parameters 
considerably. We found the effects of instrumental resolution and 
its variation along the dispersion axis to be very important. There- 
fore, we measured the resolution of our spectra in the middle of the 
available spectral range using arcs. In spite of a careful reduction 
and analysis, Figure 16 shows systematic shifts between different 
resolution data. We conclude that it is important in the future to ob- 
tain higher resolution data and to measure and model the variations 
of spectral resolution for each spectrum in the entire spectral range. 



Hot subdwarfs in the GALEX survey 17 




100000 10000 2000 



Figure 15. T c n - log(L/LE,jj) distribution of stars in the MILES library 
(light grey), WD sample (logg > 6) from SDSS DR4 (Eisenstein et al. 
2006, dark grey) and the GALEX subdwarf sample (black). The EHB, 
HeMS and WD sequences and three luminosity classes from Kurucz (1993) 
are marked in the figure. Point sizes are proportional to the inverse of sur- 
face gravities (point size = 2/(logg+ 1)) for the sake of clarity. Composites 
from the GALEX sample are marked with crosses. 
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Figure 16. Systematic shifts in the temperature-gravity diagram between 
NTT EFOSC2 observations at different resolutions and/or spectral cover- 
age. Data points of the same stars are connected. We measured A/1 = 6.4 
A for grism #7 (3700 - 5200 A) and AA = 13.7 A for #1 1 (3700 - 7400 
A). The different spectral window and low-order variations in the resolution 
can result, at worst, in AT^g a 5000 K and Alogg m 0.6 systematic shifts. 
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Figure 17. T e g - logg diagram indicating different instruments and res- 
olutions for all of the 191 observations. There is no correlation between 
instrumental setups and derived atmospheric parameters for the higher res- 
olution data, but only less accurate parameters could be derived from the 
low-resolution observations. 

A combination of the x 2 and residual minimisation (or equivalent 
width fitting) will be implemented in XTGRID to overcome the 
small scale variations caused by spectral convolutions. We found 
that high series members of H lines (from H, 5 to the Balmer-jump) 
are important in temperature and gravity determination and data 
quality is essential in this spectral range. Figure 16 shows that the 
observed spectral range can also be a source of systematic devia- 
tions, especially the inclusion of the H„ and H^ lines, which form 
higher in the atmosphere and may show non-LTE effects. 

Systematic shifts that arise from different resolutions or differ- 
ent spectral coverage are represented in Figure 16 where we com- 
pare the parameters derived from spectra taken with grism #7 and 
#1 1 with the EFOSC2 spectrograph on the NTT. Shifts from the in- 
clusion of the H ff line or due to the strong blending of spectral lines 
in low-resolution spectroscopy can be seen. This problem affects 16 
observations of 14 stars in our sample, seven of these stars, shown 
in Figure 16, have higher resolution observations as well, hence we 
regard the effects on our final results negligible. The higher res- 
olution NTT and Mayall spectra provide consistent atmospheric 
parameters as shown in the temperature-gravity plane in Figure 
17. We conclude that spectral resolution affects the derived atmo- 
spheric parameters and low resolution data provides systematically 
lower temperatures and gravities for sdO stars. These trends are not 
so remarkable in sdB stars, but considerable shifts can be found. For 
detailed subdwarf modelling in the optical, not just a high signal- 
to-noise, but medium to high resolution data is necessary (A/1 < 3 
A, or R > 3000). 



7 SUBDWARF POPULATIONS 

The observed distribution of subdwarf stars can be explained in 
the framework of both the canonical formation theory (Han et al. 
2003a, Han et al. 2003b, Mengel, Norris & Gross 1976) and the 
hot-flasher scenario (Miller Bertolami et al. 2008, Lanz et al. 2004, 
D'Cruz et al. 1996), with an interplay between atmospheric diffu- 
sion and stellar winds (Unglaub 2008). Both theories have three 
formation channels that we describe in the context of our observa- 
tions. 

In Figure 1 8 we present evolutionary sketches in the r efr - He 
diagram of the canonical and hot-flasher formation models, and in 



Figure 19 we show predictions if diffusion and stellar wind may 
change the observed He abundance. 

First, we describe the branches of the canonical formation. 
The common-envelope channels (Han et al. 2003a) predict subd- 
warfs with low-mass MS (CE-1 channel) or WD (CE-2 channel) 
companions in close binaries with 0.1-10 days orbital periods. 
Such binaries are usually discovered in radial velocity surveys, be- 
cause in these systems the low-mass companion is hardly visible 
next to a bright subdwarf. These subdwarfs are predicted to have 
a very thin H-rich envelope and a sharp mass distribution peaked 
around 0.46 M . Our observations confirm a population near 28 000 
K and logg = 5.45. 

Hotter and more helium abundant sdB stars near 33500 K 
and logg = 5.8 are predicted by the Roche-lobe overflow channel 
(ROFL-1, Han et al. 2003a). These stars are expected in wide bina- 
ries with 1 -2 M G companions, 400- 1500 days orbital periods and 
with thin H envelopes. The location of possible rapid pulsators in 
our survey coincides with this population and shows a higher com- 
posite frequency. Interestingly, the rapid pulsator class-prototype, 
EC 14026-2647 itself, and the three rapid pulsators discovered later 
are all in very similar composite spectra binaries (O'Donoghue et 
al. 1997). 

He-sdO stars may form by WD binary mergers (Zhang & Jef- 
fery 2012). This channel would intrinsically explain single subd- 
warf stars that are expected to show He dominated atmospheres en- 
riched in CNO-processed nuclei. Following the He core-flash these 
stars evolve through small loops (few 1000 K) in the temperature- 
gravity diagram of Zhang & leffery (2012) and reach their locations 
at various temperatures depending on their C and N abundance. 
Binary merger channels have been proposed with low-mass stars 
as well as with substellar companions (Soker 1998). The merger 
theory assumes a large rotational velocity of the progeny that is 
not generally observed in single subdwarf stars. Geier et al. (2009) 
found v mt sin/ < 10 km s -1 for ~100 single sdB stars. However, 
Lanz et al. (2004) found large rotational velocities for two C-rich 
sdB stars (PG 1544+488 and JL 87) and Hirsch, Heber & O'Toole 
(2008) found two out seven sdO stars having a projected rotational 
velocity in excess of 20 km s -1 . Recently, Geier, Classen & Heber 
(201 lb) found EC 22081-1916 a single sdB star that is a fast rota- 
tor with v rot sin i « 150 km s -1 and most probably a merger prod- 
uct of a WD and a substellar companion. In close binaries tidal 
spin-up can cause such high rotational velocities. [Note added in 
proof: After this paper has been accepted, Vennes et al. (2012) re- 
ported that one of our targets (11411-3053) is a fast rotating sdB 
star (VrotSini = 164 ± 5 km s ) in a short period sdB-WD binary.] 
A targeted radial and rotational velocity study would be helpful to 
constrain the properties of our subdwarfs and investigate the merger 
channel. 

The hot-flasher scenario (Miller Bertolami et al. 2008) also 
predicts three major groups of subdwarfs. In the early hot-flasher 
channel the core-flash occurs after the star leaves the RGB, during 
the constant luminosity phase and hot subdwarfs are predicted with 
standard H/He envelopes. Possible rapid pulsators in our catalogue 
show such atmospheres; they are He enriched and have notably less 
CNO in their atmospheres. The late hot-flasher (shallow mixing) 
theory assumes a core-flash early after the star reached the WD 
cooling track. This concept predicts a He and N enriched H atmo- 
sphere, similar to possible slow pulsators or cooler He-sdO stars 
in our sample. The late hot-flasher (deep mixing) predicts strong 
He, C and N enrichment, exactly as we found in the case of hot- 
ter He-sdO stars. All these features can be observed in Figure 8 
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Figure 18. Possible evolutionary sketches for sdB, sdO and He-sdO stars based on our observations and the predictions of the canonical (left) and hot-flasher 
(right) scenarios. 



and 9 as well. In the hottest sdO stars we found signs of increasing 
abundances of C and N with temperature. 

Since the evolution from the primary He core- flash to the EHB 
or HeMS is about 40 times shorter than the evolution on the EHB, 
stars are expected to group around the EHB or the HeMS. Such an 
accumulation is not observed for He-sdO stars. This may suggest 
that these subdwarfs evolve through the HeMS in large numbers 
without slowing down, which is physically impossible. Hence we 
conclude that they show more diverse atmospheric parameters be- 
fore settling on the EHB as sdB stars or starting their evolution 
towards the WD cooling sequence. 

The left panel of Figure 18 shows evolutionary sketches in the 
T e g - He diagram based on our observations and the predictions 
of the aforementioned canonical studies. The distribution of sdB 
and He-weak sdO stars can be explained according to the canoni- 
cal theory. Convection and rotation on the HB are significant below 
1 1 500 K. BHB stars show a decreasing He abundance with sur- 
face temperature until core He burning is on. After core exhaustion, 
these stars evolve towards the AGB and lower surface temperatures. 
Along this evolution the He abundance is expected to increase as 
shown by lines a and b in Figure 18. In sdB stars He is sinking 
to about 24 000 K (O'Toole 2008). Over 24 000 K the increasing 
UV flux starts to bring He to the surface and steadily increases the 
He abundance until about 36 000-38 000 K where core He burning 
stops. In stars with a lower He abundance (thicker H envelopes and 
larger total mass) this happens at lower temperatures. Canonical 
models then predict a fast evolution to higher temperatures. After 
passing a maximum temperature, stars cool and He sinks again, 
these post-EHB stars rapidly evolve to WDs. We sketched and la- 
belled this evolution with lines c, d and e. 

The canonical formation of He-sdO stars (Zhang & Jeffery 
2012) predicts different channels depending on the merger type. 
Assuming He-WD binaries it predicts C and N enriched He-sdO 
stars in agreement with our observations. In Figure 8 He-sdO stars 
separate in temperature according to their C and N abundance. If 
these stars belong to the HeMS, this separation suggests a differ- 
ent average stellar mass as well. Following the He core-flash the 
slow merger model includes accretion from a debris disc and pre- 
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Figure 19. Evolutionary sketches for sdB, sdO and He-sdO stars based on 
our observations and theoretical predictions with atmospheric mixing and 
stellar winds. 



diets N-rich low-mass He-sdO stars (line g). The fast merger model 
involves a period of coronal accretion and predicts higher mass C- 
rich subdwarfs (line h). The composite model predicts a similar 
abundance pattern, but a broader mass distribution where low-mass 
He-sdO stars are enriched in N and high mass stars are enriched in 
C, hence this evolution proceeds between the slow and fast merger 
tracks. After core He exhaustion they may evolve directly to lower 
temperatures and abundances as indicated by the dashed lines ( and 
j, or reach a maximum temperature before joining the WD cooling 
tracks (line/). 

The right panel of Figure 18 shows evolutionary sketches 
based on the hot-flasher scenario of Miller Bertolami et al. (2008). 
This predicts an evolutionary link between He-sdO, the hotter sdB 
and He- weak sdO stars (line m). The different binary frequency of 
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He-sdO and sdB stars implies a mixed sdB population, therefore 
we repeat the sketches for the canonical evolution (lines c, d and 
e) that may contribute to sdB stars. The early hot-flasher and late 
hot-flasher (shallow mixing) channels may also contribute to the 
mixed sdB population. The hot-flasher formation of He-sdO stars 
predicts larger loops (dashed lines) from He shell-flashes than the 
canonical theory, extending to 10 000 K. These loops might cover 
the entire temperature range of He-sdO stars in our sample. Along 
these loops stars may cross the wind limit and show changing en- 
hancements of C and N. In light of the hot-flasher scenario, our 
observations imply that stars after a He shell-flash may show C 
enriched atmospheres (/) that may transform to N-rich CNO pro- 
cessed atmospheres in their cooler phase (k). After these transients 
stars either move towards higher temperatures (line/) or settle on 
the EHB among suspected rapid pulsators and join the canonical 
evolution of the hotter sdB stars (line c). 

Between 36 000 and 40 000 K a photospheric convective mix- 
ing is predicted over logg a 5.9 increasing the He abundance un- 
til no H is measurable (Groth, Kudritzki & Heber 1985) as shown 
by line n in Figure 19. This mixing makes the tracking of forma- 
tion theories difficult and might be responsible for the atmospheric 
properties of He-sdO stars. These stars may cross the wind limit 
periodically and transient He enhancements (dash-dotted line) may 
occur while stars drift to higher temperatures. With this evolution 
to higher temperatures the C/N ratio is steadily increasing accord- 
ing to our observations. Convection ceases around 45 000-52 000 
K and He sinks again (O' Toole 2008). As a result, He-sdO stars 
may return to the canonical evolution indicated by line c through 
dashed line p, or reach a maximum temperature and follow line /. 
The gap between 40000 and 54000 K at [He/H]. « -1.5 coin- 
cides with the He-sdO region in our sample suggesting that most 
of the He-rich sdB stars may either pass through very fast (dashed 
line q), or suffer a huge metal abundance increase upon entering 
this region. A larger sample could help to statistically confirm the 
existence of this empty region. A period of deep atmospheric mix- 
ing can also explain the high 3 He concentrations observed at some 
subdwarf stars (Heber 2009). Due to a high temperature gradient in 
subdwarf atmospheres the pp-chain can produce this isotope rela- 
tively close to the photosphere. Because the equilibrium abundance 
of 3 He is much higher at low temperatures, 3 He produced in deeper 
layers may accumulate and replace 4 He in the atmosphere if diffu- 
sion brings it to the surface. 

Progenitors of low-mass WDs evolve through the parameter 
region occupied by subdwarfs. Figure 20 shows evolutionary tracks 
for different stellar masses taken from Driebe et al. (1998). Evolu- 
tional time-scales are inversely proportional to line widths in order 
to emphasise regions where an accumulation of stars is expected. 
Such low-mass pre- WDs are conceivable in our sample, in partic- 
ular among He-weak sdO stars, but, in general, the distribution of 
our subdwarfs does not confirm a significant fraction of such low- 
mass, non-core-He burning stars. 

These interpretations can qualitatively relate our observations 
to the various subdwarf evolutional theories, but deeper investiga- 
tions are necessary on a larger sample to estimate the contribution 
of the channels of different subdwarf formation theories as well as 
the role of diffusion and mass-loss. 



8 CATALOGUE OF SUBLUMINOUS STARS 

The catalogue in Table 2 and 3 contains 191 observations of 180 
subluminous field stars sorted according to increasing right as- 
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Figure 20. Low-mass He WD evolutional models from Driebe et al. (1998). 
Tracks for 0.259, 0.300, 0.331 and 0.414 M G are overlaid and labelled with 
the stellar mass. The 0.259 M Q track crosses this region twice following a 
He shell-flash. 



cension. In resolved binaries the primary and the companion are 
labelled with A and B suffixes, respectively. Our observations 
are typesetted with bold and MILES templates with normal font. 
We list GALEX Nw and GSC2.3.2 V magnitudes taken from the 
GSC quick- V photometry. Atmospheric parameters like effective 
temperature, surface gravity, He and CNO abundances are listed 
for subdwarfs and only effective temperature, surface gravity and 
metallicity for companions. The interstellar reddening was esti- 
mated towards each target and we list the values of the derived 
E(S - V) colour excess. The V— J and J-H colour indices were cal- 
culated from 2MASS photometry. In parenthesis, we give synthetic 
colour indices calculated from TLUSTY SEDs and HILIB templates 
in a similar fashion like in Figure 13. Because of the higher un- 
certainties in the V magnitudes the J-H colour indices are more 
reliable. Our comments are abbreviated with flags and these are 
explained in Table 6. We performed a cross-correlation of our sub- 
dwarf list with SIMBAD 6 and The Subdwarf Database 7 (0stensen 
2004, Ostensen 2006). In SIMBAD a search radius of 1 arcminute 
was used. Resolved stars are marked with one of their common 
identifiers in the last column. Ambiguous identifications are ap- 
pended with a question mark. Out of the 1 80 stars we found 32 in 
the PG, 16 in the First Byurakan Survey (FBS, Mickaelian 2008), 
6 in the HS, 3 in the Edinburgh-Cape (EC, Kilkenny, O'Donoghue 
& Stobie 1991) and 2 in the Hamburg-ESO Survey (HE, Wisotzki, 
Reimers & Wamsteker 1991) catalogues. We cross-correlated our 
data with the list of subdwarf pulsators in Ostensen et al. (2010) as 
well, but no match was found. A cross-correlation with Geier et al. 
(2011a) in search for close-binaries revealed only the two known 
short-period binaries: J032 1+4727 and J2349+3844, already dis- 
cussed in Kawka et al. (2010). 



8.1 Comments on individual stars and overlaps with other 
catalogues 

The new analysis confirms our previous results on PG 0314+146 
(He-sdO), PG 0838+133 (He-sdO), PG 1432+004 (sdB) and PG 



http : //simbad . u- strasbg . fr/simbad/ 
7 http : //catserver . ing . iac . es/sddb/ 
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Table 6. Flags that are used in the comments of Tables 2 and 3. 



Flag Meaning 



c - Contamination, scattered light from a nearby star, 

t - See text; further details can be found in Section 8.1. 

#7, #11 - Spectral resolution; 

#7: ESO/NTT/EFOSC2, grism #7, AA = 6.4 A. 

#11: ESO/NTT/EFOSC2, grism #11, AA = 13.7 A. 

If resolution is not given: NOAO/Mayall/RC, grating 

KPC-10A, A/1 = 4.95 A. 
d - Binary decomposition degeneracy: components are 

either similar type, or a high metallicity single line 

binary. 

b - Spectral signatures of a cool companion. 

rv - Radial velocity with respect to the kinematic LSR is 

over 100 km s _I . 
RV - Radial velocity curve proves a close binary, references 

are given in footnote, 
f - Poor fluxing below 4500 A. 
ca, na - Can H&K or Nai D lines; indicating either a cool 
companion or interstellar absorption, 
he - He i/He n line inconsistency. 
M - Data is taken from the MILES library. 
C&N - C&N A46A0 - i4665 A blend, 
s - Possible slow pulsator 

r - Possible rapid pulsator 



1629+081 (sdB). The new fitting method provides consistent re- 
sults with Paper I. 

J0047+0958 (HD 4539) was observed with NTT and is one 
of the two subdwarfs in the MILES library observed with the 2.5m 
Isaac Newton Telescope at 2.3 A resolution. It is the brightest star 
in our sample with V = 10.24 magnitude. Cenarro et al. (2007) 
determined r cff = 25 200 K and logg = 5.40 in good agreement 
with our r eff = 24650+^° K and log g = 5.38 + _°f y Then, we refitted 
the MILES spectrum with our method and found 7/ cff = 24 000+™ 
K and logg = 5. 21+jj^. From the higher resolution data we derived 
slightly lower temperature and gravity for this star. 

J0059+1544 (PHL 932) was suspected to be the central star 
of a planetary nebula. Therefore, it has been targeted by numerous 
studies to reveal its connection. Frew et al. (2010) derived = 
33 490 ±73 K, logg = 5.81 ±0.02 and [He/H]. = -1.58 ±0.03 and 
concluded it is an sdB star exciting a H n region in the dense ISM. 
Napiwotzki (1999) derived r eff = 35 000 ± 900 K, logg = 5.93 ± 
0.12 and [He/H]. = -1.53 ± 0.05 and a mass of 0.28 ± 0.01M G . 
Our analysis confirms a slightly lower He abundance, we measured 
T cff = 33530!™ K, logg = 5.83!°^ and [He/H]. = -1.69!°;°*. 

J0321+4727 and J2349+3844 have been analysed before by 
Kawkaetal. (2010). For J032 1+4727 they obtained T cff = 29 200± 
300 K, logg = 5.5 ± 0.1 and [He/H]. = -2.6 ± 0.1 in agreement 
with our measurements: T cS = 27 990^ K, logg = 5.34 ± 0.07 
and [He/H], = -2.52^^, although our new temperature is slightly 
lower. For J2349+3844, we derived a lower temperature. Fitting 
only H„, Kawka et al. (2010) measured T cff = 28400 ± 400 K, 
logg = 5.4 ± 0.3 and [He/H], = -3.2 ± 0.1, while we measured 
r eff = 23 770!g» K, logg = 5.38!";^ and [He/H]. = -3.44^. 
Both stars show radial velocity variations and are confirmed bina- 
ries. Phase variation is a possible source of discrepancies in mea- 
sured atmospheric parameters. 

J0507-2802 (HE 0505-2806) and J0657-7324 (CPD-73 420). 
For both stars we have multiple observations showing variable 
composite spectra and a significant IR excess. Changing compos- 



ite features are also visible in their spectra. An inspection of their 
fields revealed crowding; both stars are visual binaries. The com- 
posite spectra is most probably the result of contamination from the 
nearby, but possibly independent star. Atmospheric turbulence can 
scatter light in the slit from these stars depending on observing con- 
ditions. In such cases the decomposition is important to remove the 
contamination and improve the subdwarf parameter determination, 
although, it is difficult due to the variability and chromatic aberra- 
tion of the scattered light. Reliable parameters for the nearby star 
cannot be derived. 

J0639+5156 was observed by Vuckovic et al. (2012) in a pho- 
tometric campaign to look for bright pulsating subdwarfs. This star 
is the brightest rapid pulsator to date with a main pulsating period 
of 260 s and amplitude of ~ 13.5 mmag in the R filter. A slow pul- 
sation mode with ~2 mmag at 1765 s period suggests the star is a 
hybrid pulsator. 

J0806+1527 was observed by Baran et al. (2011) and found 
to be a hybrid pulsator with two rapid and two slow modes. 

J0851-1712 (TYC 6017-419-1). We derived 7/ cff = 39 060 K 
and log g = 5.79 for the star, placing it at the end of the EHB on the 
HeMS. We did not detect H in its atmosphere while the C/N ratio 
remained below 0. 12 and the O/N ratio below 0.54 by number with 
increasing He abundance. We classify this star as a N-class object. 
The C/N ratio of 0.05 and O/N ratio of 0.2 by number is typical 
for CNO processed material. Peculiar abundances can be the re- 
sult of mixing, diffusion and mass-loss processes (Unglaub & Bues 
2001). The atmospheric parameters of the star correspond to the re- 
gion of photospheric convective mixing in the - log g diagram. 
High resolution and high signal-to-noise observations could help 
improving the CNO abundance determination of this exotic star. 

J0934-2512. Our previous analysis in Paper I found peculiar 
He i /14471/Hen /14686 line ratios inconsistent with other Balmer 
line spectra with similar He abundances. Hen lines are relatively 
strong for the temperature of the star. We verified that the rele- 
vant atomic data (energy levels, oscillator strengths, line widths) 
in our model atoms and line lists are up-to-date for these He lines. 
Therefore, the cause for the observed line strength inconsistencies 
remains to be determined. Our new analysis showed that the star is 
in the post-EHB evolutional phase in the r eff - log g diagram, but 
belongs to the possible short period sdB pulsators in the T c g - He 
diagram. The spectrum indicates a radial velocity that does not ex- 
clude binarity. To test for a subdwarf companion we repeated our 
fit starting with a single star model and included a companion from 
the spectral library compiled from all other subdwarfs modelled in 
this work. We found the spectrum might be dominated by an sdB 
star and the He n line can be the signature of a hot, but faint com- 
panion. An underestimated temperature, gravity and He abundance 
could have similar effects, but would affect the entire spectrum. The 
true nature of the He line and surface gravity inconsistency would 
require high resolution spectra and radial velocity measurements. 
To a lesser extent such inconsistencies can be observed in six other 
sdB stars with various atmospheric parameters. Because of a possi- 
ble degeneracy in the spectral decomposition of such binaries they 
are listed as single stars with flag "he" in Table 2 referring to He 
line inconsistency. 



Table 2. Atmospheric parameters of hot subdwarf stars in the GALEX survey. 



to 
to 



GALEXJ 



Type 



Mjv 
(mag) 



V 

(mag) 



T c ff 
(K) 



logs 
(cm s~ 



Abundances, [X/H], or [Fe/H] 



E(B - V) V -J J -H Comments Other name 



He 



003531.98+245917.22 


sdB 


13 


.61 


14 


16 


003639.99+375553.37 


sdB 


12 


.75 


14 


.07 


004729.22+095855.69 


sdB 


10 


.44 


10 


.24 


004729.22+095855.69 












004759.62+033742.88A 


sdB 


1 1 


.86 


12 


,39 


004759.62+033742.88B 


F6V 










004917.26+205640.85 


sdB 


13 


.68 


14 


.72 


005956.67+154413.75 


sdB 


1 1 


.95 


12 


.27 


011103.81+194743.73 


sdO 


13 


,35 


14 


46 


011525.92+192249.86A 


sdB 


13 


10 


13 


.09 


01 1525. 92+192249. 86B 


F2V 










011554.22+261401.11 


He-sdO 


13 


.84 


14 


,51 


011627.22+060314.22A 


sdB 


13 


.58 


13 


.24 


011627.22+060314.22B 


F6V 










011928.88+490109.39 


He-sdO 


12 


,63 


13 


,56 


014307.50 383316.30 


He-sdO 


1 1 


.79 


13 


.00 


015421.98+490521.91 


sdB 


13 


10 


13 


.51 


020447.14+272903.69A 


sdB 


13 


.41 


13 


.84 


020447. 14+272903.69B 


GOIV 










020656.08+143857.08 


sdB 


13 


,79 


13 


.41 


021021.86+083058.90A 


sdB 


13 


,80 


13 


,40 


021021. 86+083058.90B 


F2IV 










021618.94+275900.81 


sdB 


13 


,79 


13 


,55 


022454.87+010938.81A 


sdB 


12 


.27 


12 


.32 


022454.87+010938.81B 


F4V 










022512.44+234822.36 


sdO 


13 


.02 


13 


.94 


023251.96+441124.96 


sdB 


13 


.74 


14 


,30 


025023.70 040611.00 


sdB 


13 


.43 


13 


.02 


031738.02+144624.00 


He-sdO 


12 


,97 


12 


67 


031738.02+144624.00 












032139.63+472718.83 


sdB 


1 1 


.91 


11 


.70 


040105.31 322346.01 


sdB 


11 


.03 


11 


,20 


040105.31-322346.01 












042034.80+012041.00 


He-sdO 


12 


,51 


12 


.52 


050018.90+091204.60 


sdB 


13 


.67 


14 


.63 


050720.16-280224.80 


sdB 


11 


.94 


12 


,39 


050720.16-280224.80 












050720.16-280224.80 













C 


N 


O 




(mag) 


(mag) 


(mag) 






-4.94 > 


-4.45 > 


-4.19 


> 


0.02 


-0.474 


-0.152 




PG 0032+247 


-5.10 > 


-4.00 > 


-4.79 


> 


0.01 


-0.709 


-0.097 




FBS 0033+376 


3 87+0.33 

-3 89+ 013 

-0.22 

-4.69 > 


3 as+0-27 

-3 67+ - 14 

J -°'-0.26 

-3.61 > 


-3.81 
-3.92 
-4.13 


> 
> 
> 


0.01 
0.00 
0.00 


-0.576 
0.510 


-0.123 
0.251 


#7,t 
M.t 

#7 


HD 4539 
HD4539 
PB 6168 


[Fe/H] = 


-0.55 






0.00 


(0.583) 


(0.244) 


M 




-4.16 > 


-4 05+°- 35 
H - u -0.58 


-4.11 


> 


0.00 


-0.371 


-0.062 


#7,s 


PG 0046+207 


-4.56 > 


-3.96 > 


-4.73 


> 


0.00 


-0.426 


-0.122 


#7,t,r 


PHL 932 


-4.77 > 


-3.51 > 


-4.28 


> 


0.00 


-0.846 


-0.193 


#7 


PG 0108+195 


-4.02 > 


-3.68 > 


-3.48 


> 


0.03 


0.428 


0.086 


#7 


CHSS 3354 


[Fe/H] = 


-0.39 








(0.364) 


(0.000) 


M 




-3 10 +0M 
J • -0.79 


-2.44 > 


-2.15 


> 


0.10 


-0.468 


-0.082 




PG 01 13+259 


-3.55 > 




-3.74 


> 


0.06 


0.926 


0.169 


#7 


PB 6355 


[Fe/H] = 


-0.15 








(0.735) 


(0.177) 


M 




-2.95 > 


-2 76+ - 92 

^■'"-0.21 


-2.36 


> 


0.04 


-0.244 


-0.257 


C&N 




-2 07+ 125 
-4.67 > 


9 c-7+1.19 
— -0.59 

-4.20 > 


-1.87 
-4.86 


> 
> 


0.03 
0.07 


-0.771 
-0.332 


-0.074 
0.096 


#7,C&N 


SB 705 

LSPM JO 154+4905 


-3.15 > 


-3.20 > 


-3.25 


> 


0.05 


0.335 


0.229 


d 




[Fe/H] = 


-1.12 








(0.137) 


(0.233) 


M 




-4.04 > 


-4 47+0-33 
H - -0.61 


-4.59 


> 


0.04 


-0.464 


-0.102 


#7,ca 


CHSS 3497 


-3.74 > 


-4.15 > 


-3.50 


> 


0.03 


0.566 


0.150 


#7,r 




[Fe/H] = 


-0.34 








(0.610) 


(0.145) 


M 




-3.88 > 


_3 no+0.35 


-3.33 


> 


0.03 


-0.321 


-0.127 


na 




-3.24 > 


-3.74 > 


-3.59 


> 


0.05 


0.684 


0.173 


#7,s 




[Fe/H] = 


-0.53 








(0.564) 


(0.165) 


M 




-4.72 > 


-4.03 > 


-3.42 


> 


0.13 


-0.429 


-0.092 


#7,ca 


HS 0222+2334 


-4.94 > 


-4.19 > 


-4.37 


> 


0.06 


-0.555 


-0.108 


na.r.rv 


FBS 0229+439 


-3.49 > 


-3 89+ 038 

-0.73 


-3.68 


> 


0.17 


-0.824 


0.268 


#7 


HE 0247-0418 


9 9C+I.I8 
■ -0.20 
_ 9 5 o+1.19 
-0.53 

-3.94 > 


9 :>s+0.86 

_ 9 5J+1.18 
-0.68 

-4.30 > 


-1.54 
-1.62 
-3.88 


> 
> 
> 


0.15 
0.20 
0.00 


-0.170 
-0.107 


-0.074 
-0.052 


#7,C&N 
#11, C&N 
RV°,s,t 


PG 0314+146 
PG 0314+146 
TYC 3315-1807-1 


a 1 -7+0 27 

~ 4A > oil 

-3.53 > 


-> n~) +0 2 1 

J -"-0.27 

-3.07 > 


-4.62 
-3.10 


> 
> 


0.05 
0.00 


-0.594 


-0.143 


#7 

#11, rv 


CD-32 1567 
CD-32 1567 


_1 5Q+0.93 
'■•"-0.14 

-3.96 > 


9 9 i +0.83 
^ -0.76 

-3.40 > 


-1.80 
-4.06 


> 
> 


0.09 
0.17 


-0.435 
-0.083 


-0.140 
-0.238 


#7,C&N 
#7,rv 


HS 0457+0907 


-4.54 > 




-4.46 


> 


0.09 


1.304 


0.214 


#7,ca,t 


HE 0505-2806 


-2.87 > 


-2.43 > 


-2.07 


> 


0.00 






#ll,ca,c,t 


HE 0505-2806 


-2.97 > 


-3.12 > 


-2.47 


> 


0.15 






#ll,ca,c,t 


HE 0505-2806 



PS 

r 

s 

S 



6140 
27520+ joo 

33530! ™ 

42780+ gjj 

31060+ 680 
jiuuu 65() 

6710 

5441 n+ 4750 
35310+ 620 
6170 



,+ 740 
- 700 


5 


.66+0; 


,+ 500 


5 


.89+0; 


-1330 




,+ 590 

- 200 


5 


.38+0 


,+ 180 

- 380 


5 


.21+0; 


,+2250 

- 970 


6 


.14+0; 



32100; 
29260; 
5500 



6580 

25480+ «o 
28000:»« l 
6430 

46670+^0 
33260+ «J 

28560! 'to 
47520+1030 

H/ J - 820 

45310+ 1610 

,JJ1 -1410 
27990+ 460 

30490! |g 

27050+"'° 

- 910 

- 510 

36270+ 490 
-1130 
180 

- 200 
25710+3450 

26070 



46050+ 



29840+ 



-1450 
+2430 
-1040 



4.14 

5 55+0- 07 

J -0.06 

= 83+0.04 
J -0.05 
c Q7+0.09 
J -°'-0.07 

6 12+0- 20 

-o.i6 

4.38 
5.89! 



,+0.17 
-0.24 
6.08!0; 24 



4.09 



,+ 510 


5 


Sfi+007 


- 500 


■ 86 -0.21 


,+«)70 


5 


64+0.06 
-0.58 


- 140 


,+ 410 

- 240 


5 


50+0.32 
J -0.02 


,+ 1320 
- 630 


5 


6 2+o.i4 

oz -0.16 



3.08 



,+ 660 

- 80 


5 


77+0.05 
' -0.06 


,+1380 
- 990 


5 


74+0.27 
-0.17 



5.45+ 



4.19 

c o it +0.04 
-0.05 

r+0.18 

-0.10 

4.00 

c a-!+o.io 
5.83_ 012 

c 7^+0.09 
■ -0.10 



5.67+ 



7+0.11 
-0.09 

5 96+0 .15 

J - 7U -0.19 

5 74+0- 27 

J -' -0.47 
5 34+0.07 
J.J+„ 07 

c 7 1 +0.06 
J - -0.04 
5 35+ 017 

c 07+O.I6 
J - y '-Q.U 

5 75+ 015 

■ -0.13 
I +0.03 
-0.06 
c 98+0.40 
•'■ 28 -0.26 

5 36+0- 22 

—0.31 



5.61 + 



-2.91 > 

-3.20 > 
_ 2 49+0.20 

-0.07 
_9 9^+0.07 
-0.09 

-2 63+0- 44 



+0.16 
0.23 

3+0.06 
-0.04 
9 97+0.17 
-z.z/„ 17 

-2.57 > 



-2.48 
-1.69 



0.10+ 



-,+0.60 
J -0.06 

-1 70+0 13 

'•'"-0.13 



16+o 44 

ulu -0.05 



1.00 



-1. 
-2.34+ 



0.36 
+0.15 
-0.08 



-2 61+ 015 

iM -0.24 

_j 4J+0.21 
1 - H1 -0.26 



9 a i +0.19 
" 2 - S1 -0.23 

-2.54 > 



-2.18+ 



d+0.24 
3 -0.24 
-1 70+ 008 



-2.56+ 



0.34+ 



t,+0.25 
-0.35 

90+o.n 

v y -0.21 
1+0.76 
-0.10 
9 C9+O.I7 
J -0.22 
Q9+O.O6 
1 - ?z -0.04 

1 81+0 25 

97+0 .12 

-0.22 

-1 46+o 14 

1 - HU -0.15 
1+0.11 
-0.14 

-1.70 > 



-2.44+ 



-2.13 



+0.49 
-0.98 
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Table 2. Atmospheric parameters of hot subdwarf stars in the GALEX survey. Continued. 



GALEXJ 



Type 



A+jv V 
(mag) (mag) 



T e ff 
(K) 



logg 
(cm s -2 ) 



Abundances, [X/H]. or [Fe/H] 


E(g — y) 


V — I 

V — J 




Comments 


Other nsmc 


He 


C 


N 


O 


(mag) 


(mag) 


(mag) 






-3 05+ - 48 
Jp -0.78 


-3.95 > 


-3.29 > 


-2.87 > 


0.10 


-0.205 


0.005 


#7,rv 




9 cn+0.45 
J -0.29 


-4.49 > 


-3.80 > 


-4.06 > 


<0.22 


-0.381 


-0.097 


f,b,na 


FBS 0550+641 


~ LiS -0.08 

-2 66+ 020 


a Q9+040 
-0.63 
-4.32 > 


-3.75 > 
-4.36 > 


-3.90 > 
-4.68 > 


0.17 
0.19 


-0.408 
-0.554 


0.027 
-0.142 


na,r 

na,b,rv,r^,t 




-3.21 > 


-4.46 > 


-4.14 > 


-4.08 > 


0.08 


1.070 


0.252 


#7,ca,t 


CPD-73 420 


-5.28 > 


-4.30 > 


-3.46 > 


-2.94 > 


0.25 






#ll,c,t,rv 


CPD-73 420 


-2.07 > 


-2.85 > 


-2.70 > 


-1.72 > 


17 






#1 1 c t rv 


r"pn— 7^ 4?n 


-0 25 +0 - 54 
-7 76+ - 22 

-0 26 

-1 56+ 033 

J -0.32 


-3 09+ 045 

JpU '-0 76 

-4.07 > 
-3.66 > 


-2.75 > 
-4 4fS+0.34 

-3.55 > 


-2.60 > 
-3.91 > 
-3.16 > 


0.09 
0.03 
0.03 


-0322 
-0.675 
-0.154 


0.000 
-0.067 
0.058 


he,C&N 

na,s 

na,d 


FBS 0658+627 




[Fe/H] 


= -0.25 






(0.075) 


(0.101) 


M 




-3 16+ - 48 

—1.03 


-4.69 > 


-4.00 > 


-4.17 > 


0.02 


-0.214 


-0.165 






74+ - 75 

04+ 056 
-0 09 

-3 40+ - 41 

J -0.87 


9 n o+0 .59 
" 2 - US -0.86 
9 79+O.85 

-0 27 

-3.99 > 


-2.09 > 

-2 76 +0 96 

-0.33 

-4.07 > 


-1.47 > 
-2.12 > 
-4.05 > 


<0.25 
0.06 
0.00 


-0.476 
-0.421 
-0.814 


0.115 
-0.197 
-0.054 


t'.rv 
C&N 

s 


FBS 0709+751 
FBS 0730+617 


i cn+o.ii 

1,7J -0 17 


-6.01 > 


-4.04 > 


-3.78 > 


0.19 


0.071 


0.407 


na,r 




9 7 c+0.13 
- -0.25 


-4.04 > 


-4.07 > 


-4.19 > 


0.01 


-0.654 


-0.110 


s 




-3.24 > 


-4.40 > 


-3.74 > 


-3.58 > 


0.17 


0.295 


0.353 


#7,rv 




-0 95+ 04 

i o<: +0.05 

-t.si_ 013 

_7 Aq+0-27 
-0 30 


-3 12+ 015 

-0 29 

-4 01+ 037 

u -0.51 

-4.04 > 


9 91+0.28 
->- z ' 1 _0.47 

-4 07 +0 - 35 

u -0.30 

-3.74 > 


-2.88 > 
-3.89 > 
-3.69 > 


<0.19 
0.00 
0.10 


-0.640 
-0.403 
-0.518 


-0.122 
-0.036 
-0.227 


#7 


FBS 0742+625 
Ton 287 ? 


47+° : «> 
"■^'-0.13 

^ -0.18 


9 7 i + 1.19 

-0.32 

-4.43 > 


-1.84 > 
-4.06 > 


-2.08 > 
-3.72 > 


0.08 
0.01 


-1.009 
-0.335 


-0.136 
-0.093 


#7,C&N 
na 


Ton 835 ? 


a c/i+0 49 


—2.27 > 


9 A A +0 22 

-2.44_ 202 


—2.32 > 


0.01 


-1.581 


-0.148 


#7 




-3.44 > 


-3.96 > 


-3.68 > 


-3.28 > 


0.04 


-0.437 


-0.117 


#7,rv 




-3.21 > 


-4.16 > 


-3.93 > 


-3.12 > 


0.02 






#7,rv 




-3.25 > 


-3.80 > 


-3.51 > 


-2.74 > 


0.08 






#7 




-2.67 > 


-3.70 > 


-3.74 > 


-3.86 > 


0.05 


-0.743 


0.377 


#7,s c ,t 




i 9 c +0.05 
-0.08 


-3 15+ 018 

J -0.15 


_3 9Q+0-25 
-0.37 


-3.04 > 


0.01 


-0.186 


-0.104 






i 79+0.11 
lw -0.16 


-4.23 > 


-3.62 > 


-4.12 > 


0.09 


-0.211 


0.534 


#7,r,rv 




-2.90 > 


-4.10 > 


a 7O+0.46 
-3-/8_ 088 


-3.64 > 


0.19 


-0.731 


-0.025 


#7 




i O7+0.17 
~ LS/ -0.22 


-4.52 > 


-3.59 > 


-3.24 > 


<0.26 


0.545 


0.424 


na,b 




i 90+0.10 

-0 62+ 05 

u - oz -0.08 

-2.62 > 


_A 4^+0.62 
H ^"-0.82 
-4.00 > 

-3.82 > 


_l 79+0.39 
J -' 7 -0.61 
n 4 7+0.25 
J - H/ -0.38 

-4.13 > 


-i C9+0.29 
J J -1.08 

-3.96 > 
-4.03 > 


0.10 
0.14 
0.04 


-0.433 
-0.645 
-0.464 


-0.112 
-0.161 
-0.111 


#7,s,rv 
#7 

#7,rv 


SDSSCGB 7585. 
TD1 32707 



050735.80+034815.00 


sdB 


13 


.82 


14 


.42 


055542.07+640812.60 


sdO 


13 


07 


13 


,33 


061325.26+342053.09 


sdB 


13 


,60 


13 


.63 


063952.00+515658.00 


sdB 


12 


.47 


11 


,97 


065736.93 732449.60 


sdB 


10 


.99 


11 


.90 


065736.93-732449.60 












065736.93-732449.60 












070123.11+532938.84 


He-sdO 


13 


.57 


15 


.31 


070331.53+623626.35 


sdB 


12 


.67 


13 


.10 


071029.43+233322.24A 


sdO 


13 


.87 


14 


.01 


071029.43+233322.24B 


F6V 










071053.44+562446.77 


sdB 


13 


,50 


14 


,98 


071618.19+750241.99 


He-sdO 


13 


.74 


15 


,04 


073527.00+613244.38 


He-sdO 


13 


,60 


14 


,18 


073804.54+262444.18 


sdB 


13 


.77 


14 


.42 


074148.99+552451.19 


sdB 


13 


.74 


14 


.71 


074435.08+302108.96 


sdB 


13 


,84 


14 


.64 


074617.10+061006.00 


sdB 


13 


,08 


12 


.85 


074722.24+622547.10 


sdB 


1 1 


.67 


11 


,90 


074855.82+304249.36 


sdB 


13 


11 


14 


16 


075147.00+092526.00 


sdB 


??? 


14 


,12 


075234.20+161604.00 


He-sdO 


13 


.87 


14 


.23 


075603.77+222628.82 


sdO 


13 


.12 


14 


,50 


075807.55-043203.87 


He-sdO 


11 


.81 


12 


,20 


080510.90-105834.00 


sdB 


12 


.35 


12 


.21 


080510.90 105834.00 












080510.90 105834.00 












080656.70+152718.00 


sdB 


13 


.74 


14 


.69 


080929.50+054843.82 


BHB 


13 


.48 


13 


84 


081203.80+070352.00 


sdB 


13 


.74 


14 


.78 


081233.60+160121.00 


sdB 


12 


,79 


13 


.57 


082556.78+130754.73 


sdB 


13 


.62 


14 


59 


082734.90+175358.00 


sdB 


13 


.72 


14 


,62 


082832.80+145205.00 


He-sdB 


11 


.65 


11 


.78 


083412.30+071211.00 


sdB 


13 


.73 


14 


86 



23990+ 630 

- 610 

45600+3™ 
34250! 33 <J 
30400+ ™ 
29940+ **j 
28930!™ 
25010+ 2 «« 

67230+S 
28750+ 3 ™ 
5 5 4 20+ 24300 
6310 

39970+ 980 
jyyiu 94() 

43200+ I™ 

48400+ 1410 
totuu 10g0 

29740+ 4 °° 

33470+ 43 ° 

27910+ 470 
z/?iu 28 „ 

26530+ ™ 

22550! *o 
30760! 83 ° 
30620! 49 ° 
49280! 5 || ° 
51160+ 580 
41360! 22a 
22320! 33 ° 

- 290 
u - 230 

22020! 340 
29410! g° 
21290! |o 

355 90+ 560 
JJJ - 620 

31580! 44 g 

38070! ^ 

29550+ 9 ^° 
i,JJU - 800 

37060! ™ 
29870! 48 ° 



22240+ 



5.73+ 



5.40+ 



<; 49+O.O8 
J - -0.11 
3+0.04 
■'-0.39 
c 7^+0.10 
-0.06 
c ro+0.13 
J - JO -0.04 
<: 4 s+0.07 
J>tJ -0.15 

5 10+ 013 

J ' lu -0.21 

4 85+ 036 
^- -0.26 

5 60+ 08 

-0 
9+0107 
-0.04 
^ 9 1 +0.43 
J - zl -0.14 

4.26 

5 9 1+0.08 
J - yl -0.14 

5 55+ 53 

J J -0.06 

6 05+ 017 

u - -0.21 

5 48+ - 06 

J - HO -0.06 

5 76+ - 04 

J - -0.12 
5 42+0-08 
J -0.05 
<: 7 1 +0.06 
' ' -0.07 

4 99+0.05 
^■"-0.06 

1+0.17 
-0.04 

5.7<» 
6.08+"' 9 

<: 71+O.O3 
J - -0.26 
<: oc+0.26 
-0.06 

5 fio+0.03 

1 +0.03 
-0.04 

5 65+ - 06 

J ' -0.06 
5 44+O.H 
3 -^ 4 -0.12 

9+0.03 
-0.05 

^+0.09 
-0.10 

5 56+ - 10 

J J -0.13 
1+0.14 
-0.10 
<: ^c+0.09 
JJ -0.08 
1+0.09 
-0.10 

5 56+°- 12 

—0. 10 



5.71+ 



5.71+ 



4.8 

5.75+ 



5.91+ 



5.92+ 



to 
S 

3> 



s 
I 

X 

R 

•3 



to 



b Vuckovic et al. (2012) 
c Baran et al. (201 1) 



Table 2. Atmospheric parameters of hot subdwarf stars in the GALEX survey. Continued. 



to 



GALEXJ 



Type 



Mjv 
(mag) 



V 

(mag) 



T e ff 
(K) 



logs 

(cm s~ 



Abundances, [X/H], or [Fe/H] 



E(B -V) V - J J - H Comments Other name 



He 



084024.05+113911.53 


B 


13 


.89 


13 


,19 


084143.80+130431.00 


He-sdO 


12 


.91 


13 


.66 


084535.54+194150.18 


BHB 


13 


,33 


13 


,34 


085137.75+244151.69 


sdO 


12 


.68 


14 


.35 


085649.30+170115.00 


sdB 


12 


.22 


13 


17 


090540.90+121228.00 


sdB 


13 


.73 


14 


,67 


092308.30+024208.00 


sdB 


13 


,39 


14 


48 


092440.11+305013.16 


He-sdO 


13 


.25 


14 


.73 


092856.10+061634.00 


sdB 


13 


.44 


14 


,29 


093448.20-251248.00 


sdB 


12 


,05 


13 


,03 


095256.60 371940.00 


He-sdO 


12 


.00 


12 


.69 


100752.00-292435.00 


sdO 


1 1 


,99 


12 


87 


101756.84+551632.41A 


sdB 


12 


.15 


12 


,30 


101756.84+551632.41B 


F5V 










102358.61+691145.65 


He-sdO 


13 


.61 


14 


,69 


102728.92+745947.26 


sdB 


13 


,89 


13 


.87 


104122.80+504420.00 


He-sdO 


13 


.69 


15 


.31 


104130.40+184210.00 


sdB 


12 


.27 


13 


.24 


104148.90 073031.00 


sdB 


12 


.49 


12 


14 


110055.90+105542.00 


sdB 


12 


,93 


14 


,20 


110541.40 140423.00A 


sdB 


1 1 


.82 


11 


.54 


110541.40-140423.00B 


F6V 










111422.00-242130.00 


sdB 


13 


.16 


12 


.68 


113128.89+283739.18 


B 


13 


,70 


13 


.71 


113633.58+750652.07 


sdB 


12 


.05 


12 


.99 


120624.37+570935.70 


sdB 


13 


.71 


14 


87 


123716.47+383541.22 


B 


13 


,84 


13 


,53 


125941.60+164827.00 


sdB 


13 


.85 


14 


37 


131512.40+024531.00 


sdO 


13 


.61 


15 


10 


135450.06+055937.32 


B 


13 


.87 


13 


,54 


135629.20 493403.00 


sdB 


11 


.72 


12 


,30 


135824.60+065137.00 


sdB 


13 


.68 


14 


,50 


140209.12+383713.12 


BHB 


12 


,49 


12 


.47 


140747.63+310318.06 


B 


13 


.81 


14 


,20 


140956.88+805520.23 


sdB 


12 


,94 


13 


87 


141055.80+093256.00 


sdB 


12 


.69 


14 


10 



Q 


N 


o 


(mag) 


(mag) 


(mag) 






-3.04 > 


-2.81 > 


-2 98+ 032 
-1 40 


0.00 


0.000 


-0.041 


n a 


NVSS J084025+1 13959 


7 57 + 1.39 
A D -0.45 


-1.55 > 


-1.52 > 


0.12 


-0.643 


-0.154 


#7,C&N 


PG 0838+133 


-4.01 > 


_7 qq+0.36 
3y -0.69 


-3.56 > 


0.04 


-0.328 


-0.085 




LB 400 ? 


—3.09 > 


—3.61 > 


-3.96_ 089 


0.04 


0.035 


0.120 


n a 




-4.44 > 


_ 4 oo+ - 28 
+.uu_j ]0 


-3.74 > 


0.14 


-0.256 


-0.079 


#7,s 




-3.87 > 


-3.52 > 


-3.66 > 


0.07 


-0.578 


0.090 


#7,s 


PG 0902+124 


-4.59 > 


-4 18+ 051 
^ lo -0.90 


-3.74 > 


0.21 


-0.583 


-0.180 


#7,ca,s,rv 


PG 0920+029 


-2.89 > 


_7 7Q+0.18 
^'-0.46 


-2.43 > 


0.01 


-0.572 


-0.328 




PG 0921+311 


-4.58 > 


-3 60+ 042 

J,u -0.59 


-3.27 > 


0.17 


-0.334 


-0.121 


#7,ca,s 


PG 0926+065 


-3.84 > 


-3.24 > 


-3.23 > 


0.03 


-0.682 


-0.051 


#ll,he,t 




-1.83 > 


-1.82 > 


-1.40 > 


0.05 


-1.168 


-0.147 


#11 




-5.32 > 


-2.65 > 


-4.36 > 


0.00 


-0.148 


0.131 


#11, f 




-3.44 > 


-3.74 > 


-2.43 > 


0.00 


1.078 


0.126 






[Fe/H] : 


= -0.54 






(0.507) 


(0.132) 


M 




7 £5+0.74 

-2-65_ 044 


— Z.jy > 


—L.LL > 


<U.Uj 


—U.J / z 


— U.UUj 


1,1 v 


rrio iuzu+oy4 


-3.32 > 


-4.31 > 


-3.41 > 


<0.42 


-0.071 


0.268 


f,na 


FBS 1023+753 


7 77+0.86 

z -"-o.n 
-3.98 > 


- 2 -87t 5 6 5 

t ^7+0.24 
J - -0.73 


-2.27 > 
-4.05 > 


0.00 
0.17 


-0.343 
-0.430 


-0.224 
-0.099 


C&N 
#7,ca,r 


PG 1038+510 


-3.76 > 


7 O7+0.40 
-3.8/_ 057 


-3.75 > 


0.09 


-0.062 


-0.093 


#7 




-4.40 > 


-3.53 > 


-3.93 > 


0.15 


-0.726 


-0.014 


#7,ca,r 




-5.59 > 


-1.46 > 


-1.24 > 


0.00 


0.743 


0.161 


#11, f 


EC 11031-1348 


[Fe/H] 


= 0.00 






(0.692) 


(0.179) 


M 




_ 401 +0.62 

-1.28 


-3.47 > 


-3.29 > 


0.14 


0.452 


-0.116 


#7,ca 


EC 11119-2405 


-4.40 > 


-3.66 > 


-2.89 > 


0.00 


-0.320 


-0.128 


n a 


CHSS 1935 


-4.91 > 


-4.17 > 


-4.40 > 


0.00 


-0.795 


-0.098 


he,na,ca 


FBS 1133+754 


—4.89 > 


—4.30 > 


—4.55 > 


0.09 


—0.706 


—0. 115 


he,r 


D/"* 1 707 i 57/1 
r(j IZUj+J /4 


-4.79 > 


-3.68 > 


-3.53 > 


0.02 


-0.383 


-0.072 


rv 


HZ 30, SSII 182 


-3.88 > 


-3.75 > 


-3.99 > 


0.13 


-0.467 


-0.139 


#7 


PG 1257+171 


-4.49 > 


-6.42 > 


-3.60 > 


0.25 


-0.055 


-0.071 


#7,ca 




-2.96 > 


-2.55 > 


-2.60 > 


0.00 


-0.140 


-0.059 






-4.17 > 


-4.00 > 


-4.37 > 


0.10 


0.317 


0.278 


#7,ca,r 


CD-48 8608 


-4.39 > 


_ 4 3Q+0.39 


-4.03 > 


0.00 


-0.474 


-0.017 


ca,s 


PG 1355+071 


_7 nq+0.21 
Z '^-0.14 
-4.21 > 


-3 10+ 028 
J1 -0.43 

-4.18 > 


-2.49 > 
-3.67 > 


0.00 
0.12 


-0.123 
-0.140 


-0.163 
0.134 


ca,rv 


PG 1400+389 
BPS BS 16082-0122 


-3.38 > 


-3.84 > 


-2.58 > 


<0.20 


-0.401 


-0.054 


f 


FBS 1411+812 


-4.23 > 


-3.57 > 


-3.90 > 


0.18 


-0.594 


-0.125 


#7 


PG 1408+098 



s 

S 



13270 
47660; 
21570; 
70000 



. 40 

- 90 
,+3710 

- 770 
210 

- 170 
,+3190 
-5940 

29270! lf 
28050! 640 
29100! f 4 ° 
40830! [go 
26060! 450 
34440+ wo 

44230!' 
41880.^0 

24260!!«° 
6610 

46750+ 940 

35370+ 800 
51150+ 1080 

J " -1000 

34890+ 440 
27440! g , 
31970! 3f0 

30120!5™ 



-3790 

6430 

23430! 48 ° 
19000! Ifo 
36030; " 



t- 520 
- 590 

33620! H° 

17570!™ 
36830!™! 
41210+ 1420 

14730+ 80 

200 

33070! Z 

27680! J!!! 

19970+ 100 

390 

24900+ 50 
zt?uu_ 3 o 5 o 

25380! 3 ,' 9 ° 
36580! mo 



4 14+0.02 
-0.02 
5.68! 5 4 

4 S7+0 03 

5 56+ 07 
J - J -0.10 

5 30+O.20 
J - J -0.03 

5 50+°- 11 

c 4, +0.09 
J -0.07 

5 84+°- 15 

J - -0.04 
c 4^+0.07 
J - H -0.08 

5 17+ 021 



5.57! 



6.41 + 



-0.36 
7+0.67 
-0.71 
5 05+ 048 

4.09 

5 63+ 43 
-0.06 

5 65+ 012 

J - -0.18 
I +0.04 
-0.66 

5 63+ 009 

J - -0.06 
c Q7+0.08 
3 - S/ -0.10 

5 90+° 68 

Jy -0.82 
4.03 

5 29+0.08 

4 40+0.04 

H ' HJ -0.05 

5 qo+0.06 

5 90+° 04 

J - yu -0A3 

3 96+° 02 



5.95+ 



r+0.09 

-0.11 

5 4 6 +o.i2 

J - HU -0.13 

4 18+° 03 

H ' lo -0.07 

5 74+0.07 
J -' -0.16 

5 36+° 06 

4 81 +o .03 
^■ ol -0.08 

; +0.03 
-0.06 
c ^n+0.10 
-0.17 

5 77+O.II 



4.25+ 



1 77+O.O4 
1,z '-0.11 
73+0-34 
"■' -0.28 
-1 77+0.03 

1,/3 -ai4 

_1 7Q+0. 17 
'•"-0.14 



-2.79+ 



-2.81 > 

3+0.37 
-0.76 
-2.87 > 

59+° 54 

_7 7 n+0.23 
^•'*-0.46 

-1.85 > 

53+ 150 
U - J -0.10 

-3.73 > 
-2.00 > 



45+° 56 
u - -0.03 

1 70+0.16 

1- ' J -o.28 

,+0.35 
-0.06 
-1.68+™ 9 
4+O.I6 

-0.23 

-7 00+° 15 
z - uu -0.14 



o.ir 



-2.44+ 



-1.04 



+0.89 
1.07 



-2 46+° 19 

Z -^°-0.31 

-2 30+ 013 

-0.09 
7+0.18 
-0.16 
1 qo+0.08 
-1.88_ 012 

-2 04+° 11 
z - u ^-0.18 

£,+0.10 

J -0.13 
_ 2 7Q+0.40 
Z -' y -0.67 

-1 09+° 15 

1 - u7 -0.06 
_7 7<:+0.25 
Z - /;l -0.43 



-2.22+ 



-1.66+ 



-2.95+ 



^+0.26 
-0.64 

-0 72+° 09 
u - /z -0.04 
1 1 H+0.08 

" L1S -0.09 

-2 76+ 028 

z - ,D _0.58 
5+0.12 



-1.78 



-0.18 



Table 2. Atmospheric parameters of hot subdwarf stars in the GALEX survey. Continued. 



GALEXJ 



Type 



TVuv V 
(mag) (mag) 



(K) 



(cm s -2 ) 



Abundances, [X/H], or [Fe/H] 



He 



N 



E(B - V) 
(mag) 



V- J 

(mag) 



J-H 

(mag) 



Comments Other name 



141115 


.90 305307.00 


sdB 


1 1 


,63 


1 1 


90 


30150 


141115 


,90 305307.00 












25080 


141133 


45+703736.95A 


sdB 


12 


,71 


13 


17 


21170 


141133 


45+703736.95B 


F0.5V 










7160 


142747, 


28 2701 08 00 

had — / vlVUiVv 


sdB 


1 1 


,33 


1 1 


48 


31880 


143519 


,80+001350.00 


sdB 


12 


,25 


12 


36 


23090; 


145928 


,50+190350.00 


sdB 


13 


,33 


14 


17 


37350; 


150109 


,07+412137.20 


B 


13 


.88 


13 


55 


15530 


150346 


,81+422125.12 


sdB 


13 


,68 


14 


80 


29170 


151325 


75+645407 12A 


sdB 


1 3 


.60 


13 


95 


34310 


151325 


75+645407.12B 


GOV 










5890 


152513 


08+605321.53A 


sdB 


12 


.13 


12 


88 


34380; 


152513 


,08+60532 1.53B 


GOV 










5820 


154509 


,52+595505.20 


sdB 


13 


,62 


14 


56 


29060 


160131 


,30+044027.00 


He-sdB 


13 


,48 


14 


56 


40330 


160209 


18+072509.84A 


sdO 


13 


.63 


14 


63 


67570 




1 S+D79S0Q 84R 


GOV 










5980 


161327 


,66+085355.54 


sdB 


13 


,88 


14 


44 


31840; 


161902 


78+483144.11A 


sdB 


13 


,35 


13 


36 


34670; 


161902 


,78+483144.1 IB 


F0V 










6700 


162431 


88+40591 3 03 


sdB 


1 3 


,64 


14 


32 


27590 


162716 


,50+402728.48 


sdB 


12 


.29 


12 


60 


23560 


162822 


,57+272919.71 


sdB 


13 


11 


13 


72 


25380; 


163201 


,40+075940.00 


sdB 


13 


,41 


12 


97 


38110 


163201 


,40+075940.00 












35780; 


163215 


,73+851359.16 


sdB 


12 


,58 


12 


45 


33590; 


163705 


,28+411541.02 


sdB 


12 


.90 


14 


15 


26850; 


164614 


,55+062300.71 


sdB 


13 


,72 


14 


84 


33170; 


165054 


,44+353341.65 


sdB 


13 


,85 


15 


49 


33020 


165157 


,16+262551.56 


sdO 


13 


,56 


14 


26 


67700 


170211 


,83+483024.66 


sdB 


13 


,22 


14 


30 


26010 


171708 


,53+675712.24 


WD 


13 


,76 


13 


73 


14910 


172748, 


,43+164456.52 


sdB 


12 


.98 


13 


35 


28210; 


172704, 


,66+294229.22 


sdB 


13 


,78 


15 


03 


37610 



- 310 
,+ 1260 

- 730 
1500 
1110 

360 
290 
780 

- 250 

^ 460 
-1090 
h 80 

- 210 

- 380 

- 250 
-I6S0 
-1560 

hlOOO 

- 690 

- 280 

- 250 

- 860 

- 770 
-15540 
-3950 

- 310 

- 320 

,+2370 
1610 

430 
260 
410 

- 250 

' +3 6? 
■ 570 
680 
1420 
1580 

210 
810 

,- 420 

- 310 
240 
390 
610 

- 280 

,+2370 
-7080 
' 1110 
140 
100 

- 150 
550 
420 
990 
1010 



5.22+ 



c 79+0.04 
J ■ -0.07 
1+0.19 
-0.11 
5 55 +031 

4.15 

5 70 +0 05 

-0.08 
c 70+O.O8 
-0.08 

S Q7+0.10 
J -' -0.15 
j+0.02 
-0.05 

5 52 +012 
J - -0.04 

6 17 +030 

-0.36 

4.37 

7+0.17 



4.18 



5.77 
3.92 

5.43 + 



-0.15 



? +0.05 
-0.05 

6.16+ 018 



-0.10 

+0.24 
-0.13 



5.74 
4.05 



5.87+ 



J - -0.24 

4.44 

« 2Q+0.06 
J -0.05 
S 44+0-03 
J -0.07 
c 1 1+0.05 
J - -0.12 
c -30+O.O6 
J - JO -0.09 
5 14+012 
J - 1 ^-0.40 

5 47+0- 15 

c 47+O.O6 
J - HZ -0.05 

5 84 +0 0S 

J -° -0.07 
c 07+O.O9 

c on+006 
5.8U_ 030 

8+0.09 
-0.03 

5 67 +0 04 

J - -0.03 
7+0.09 
' -0.08 

c C9+O.IO 

J ' J -0.13 



5.18 



5.27+ 



1 co +0.03 
-1.58 Q0J 

-1.39 + °J° 

-0.14 


-4.18 > 
-3.59 > 


-4 40+ - 29 

+- w -0.67 

-3.16 > 


-4.52 > 
-2.76 > 


0.04 
0.02 


-0.986 


-0.046 


#7,ca,rv 
#ll,ca,rv 


CD-30 11223 
CD-30 11223 


-2.36 > 


-3.17 > 


-3.02 > 


-1.53 > 


0.00 


1.057 


0.109 








[Fe/H] = 


0.25 






(0.333) 


(0.069) 


M 




1 71+0.05 
-011 

7 3Q+0.18 
^■ J? -0.20 


-3.87 > 
-4 66 +0 72 

H - uu -0.77 


-4 16 +033 
-0.73 

-3 26 +0 - 36 
J - -1.02 


-4.48 > 
-2.84 > 


0.08 
0.20 


-1.049 
-0.884 


-0.140 
-0.072 


#7,r 
#7,ca 


EC 14248-2647 
PG 1432+004 


_] 41+O.O8 

-0.17 
-2.53 > 


-3.70 > 
-4.09 > 


3 77+0.28 
- -0.88 

-2.89 > 


-3.76 > 
-3.34 > 


0.25 
0.03 


-0.751 
-0.139 


-0.077 
-0.109 


#7 
rv 


PG 1457+193 
PG 1459+416 


7 oo+0. 17 
— Z.00 n ~ n 
-0.30 

1 O1+0.74 
_1 - 81 -0.30 


-3.94 > 
-2.98 > 


-0.96 

-3.12 > 


-4.52 > 
-2.82 > 


0.00 
0.00 


-0.480 
0.585 


0.085 
0.186 


he,s 


PG 1501+426 




[Fe/H] = 


-1.79 






(0.681) 


(0.307) 


M 




-1 56 +0 18 

1 -' D -0.18 


-3.55 > 


-3.46 > 


-3.55 > 


0.00 


0.597 


0.192 


r 


PG 1524+611 


[Fe/H] = 


-0.74 






(0.486) 


(0.195) 


M 




-3.44 + °- 2 « 

-0.60 

-0 53 +0 04 
"""-0.21 

-1.56+02? 

-0.35 


-4.18 > 
7 7C+0.15 

-0.38 

-2.95 > 


_4.49 + °S 

-0.90 
7 71+0.19 
z - -0.41 

-2.66 > 


-4.14 > 
-2.92 > 
-3.09 > 


0.00 
0.24 
0.04 


-0.521 
-0.398 
0.492 


-0.125 
-0.127 
0.296 


s 

#7 
#7,d 


PG 1544+601 
PG 1559+048 
PG 1559+076 


[Fe/H] = 


0.07 






(0.300) 


(0.172) 


M 




7 77+O.I5 
z - z '-0.18 


-4.28 > 


-3.97 > 


-4.40 > 


0.00 


-0.857 


0.003 


#7,ca,r,rv 


PG 1611+090 


-6.96 > 


-5.86 > 


-5.98 > 


-7.13 > 


0.00 


0.705 


0.124 


r 


BPS BS 16029-0090 




[Fe/H] = 


-0.01 






(0.494) 


(0.007) 


M 




7 co+0.14 
■ -0.30 


_3 07+0.25 
J -"-0.60 
-4.35 > 


_ 4 47+0.37 
+■ -1.11 
-3.56 > 


-3.64 > 
-3.75 > 


0.03 
0.00 


-0.323 
-0.467 


0.052 
-0.157 


rv 


KUV 16229+4106 
FBS 1625+407 


_7 4<:+0.12 
z - HJ -0.30 
7 71 +0.27 
-0.29 


-3.66 > 
-5.03 > 


-> 73+0,31 
-■' J -0.45 

-4.43 > 


-2.98 > 
-4.21 > 


0.00 
0.03 


-0.455 
0.134 


-0.137 
0.225 


rv 
#7 


PG 1629+081 


-2.15 > 


-4.20 > 


-3.72 > 


-3.89 > 


0.10 






#11, he 


PG 1629+081 


-3.38 > 


-4.53 > 


-4.31 > 


-4.52 > 


0.06 


-0.401 


-0.143 


na 


FBS 1642+853 


7 7C+0.I6 
z "' J -0.24 
_ 2 43+0.23 
Z -^ J -0.12 


-4.30 > 
-3.70 > 


a on+0.28 
^-'"-0.98 

-4.08 > 


-4.01 > 
-4.87 > 


0.00 
0.01 


-0.420 
-0.617 


-0.124 
0.100 


s,rv 

#7,ca,r,rv 


FBS 1635+413 
PG 1643+063 


1 ^2+0.03 
-0.08 


-3.81 > 


_4 17+0.31 

-0.41 


-4.15 > 


0.00 


-0.241 


-0.057 


r 


PG 1649+356 


1 go+0.31 
" L83 -0.12 


-3.28 > 


-3.74 > 


-3.56 > 


0.00 


-0.608 


-0.252 


na 




7 co+0.09 


-4.00 > 


-3.90 > 


-3.67 > 


0.00 


-0.530 


-0.155 




PG 1700+486 


-2.26 > 


-3.18 > 


-2.81 > 


-2.88 > 


0.00 


0.121 


-0.083 


RV e ,ca 




7 O7+0.30 
" 2 - 87 -0.61 


-3.67 > 


-4.11 > 


-3.94 > 


0.14 


0.024 


0.265 


#7,ca 




-2.66 > 


-4.29 > 


-3.75 > 


-3.69 > 


<0.07 


-0.650 


-0.028 


f 


LB 334 ? 



S 
I 

3' 

PS 

S 
I 

X 

R 

•3 



Vermes et al. (2011b) 



to 



Table 2. Atmospheric parameters of hot subdwarf stars in the GALEX survey. Continued. 



to 

ON 



GALEXJ 



Type 



Wuv V 
(mag) (mag) 



T c ff 
(K) 



(cm s -2 ) 



Abundances, [X/H], or [Fe/H] 



E(B-V) V-J J-H Comments 



Other name 



He 





XT 
IN 


( \ 
\ J 


(mag) 


(mag) 


(mag) 






-3.84 > 


-3.81 > 


-3.60 > 


0.16 


-0.710 


-0.062 


#7,he,s 




-3.82 > 


_4 99+O.38 
-0.57 


-4.33 > 


0.08 


-0.125 


-0.053 


#7,ca,s 




-2.21 > 


-1.65 > 


-1.63 > 


0.11 


0.591 


0.214 


#11 


TYC 2084-448-1 


[Fe/H] 


= 0.17 






(0.640) 


(0.176) 


M 




-4.13 > 


-3.58 > 


-3.20 > 


0.00 


-0.138 


-0.237 


na 




9 Q9+0. 14 
M -0.55 


-3.40!"; 32 


-3.63 > 


0.00 


-0.609 


-0.252 


r,rv 


FBS 1749+373 


-2.33 > 


~ 2 - 86 -a% 


-2.67 > 


0.07 


-0.127 


-0.513 


rv 


FBS 1753+692 


-3.58 > 


-4.70 > 


-3.98 > 


0.08 


0.766 


0.253 


#7 




[Fe/H] : 


= -0.32 






(0.722) 


(0.202) 


M 




-2.06 > 
-3.50 > 


-1 46+ 034 

-4 16+ 052 
-1.76 


9 ,7+0.57 
-0.65 

-4.11 > 


0.00 
0.17 


-0.626 
-0.412 


-0.161 
-0.049 


rv 

#7,ca,r 




-4.68 > 


-3.74 > 


-3.25 > 


0.08 


-0.329 


-0.051 


#7,ca,rv 




-3.97 > 


-3.88 > 


-3.47 > 


0.12 


0.435 


0.330 


#7,ca,t 




-0.82 > 


-0.12 > 


-0.66 > 


0.12 


-0.422 


-0.035 


#11, t 




-2.58 > 


-9 4^+0.37 
-0.25 


-2.91 > 


0.07 


-0.573 


0.040 


#7.t 


CD-51 11879, LSE 263 


-3.20 > 


-3.49 > 


-3.30 > 


0.01 


-0.494 


-0.124 


#ll,rv 


BPS CS 22936-0293 


3 37+0.70 
• -0.60 


-2.69 > 


-1.82 > 


0.15 


-0.544 


-0.082 


#7,C&N,rv 




—2.32 > 


—2.20 > 


— 1.76 > 


0.15 






ff 1 1,C&JN,1V 




-4.06 > 


-3.76 > 


-4.25 > 


0.00 


-0.384 


-0.141 


r 


FBS 1930+872 


-3 67+ a09 

J -°'-0.46 

-4.32 > 
Q4+0.43 
3y -0.65 
-4.15 > 


_4 in+0.27 
+ 1V -0.69 

-4 49+0.36 
-0.40 

_o Q3+0.34 
J - -0.56 

_ 4 4 ,+0.39 
+■+'-0.98 


-4.55 > 
-4.50 > 
-3.95 > 
-4.71 > 


0.06 
0.04 
0.11 
0.04 


-0.327 
-0.105 
-0.348 
-0.543 


0.034 
0.274 
-0.195 
-0.206 


#7,ca 
#7,ca,s 
#7,ca 
#7,ca,r 




-4.37 > 


-4.34 > 


-3.57 > 


0.00 


0.357 


0.237 


#7,rv 




[Fe/H] : 
_ 4 <; 6 +0.61 
+■ -0.08 
[Fe/H] 


= -1.11 
+' D1 -0.08 

= 0.13 


-4 25+ 052 
+■ -0.08 


0.02 


(0.558) 
0.702 
(0.495) 


(0.241) 
0.266 
(0.196) 


M 

#7,rv 
M 




-5.93 > 


-3.77 > 


-3.52 > 


0.07 


0.847 


0.240 


#7 




[Fe/H] : 


= -0.25 






(0.629) 


(0.175) 


M 




-4 16+ - 29 
+' lu -0.46 


_4 19+0-29 


-4.35 > 


0.05 


-0.947 


-0.095 


#7,ca,rv 




-3.22 > 


-2.94 > 


-2.62 > 


0.13 


1.502 


0.454 


#7,rv 




[Fe/H] : 


= -0.29 






(1.289) 


(0.432) 


M 




-4.34 > 


-4.18 > 


-4.40 > 


0.01 


-0.945 


-0.111 


#7,rv 




9 -17+O.85 
"-0.22 
-4.74 > 


_9 4^+0.56 
z - -1.02 
-3.93 > 


-1.92 > 
-4.37 > 


0.03 
0.07 


-0.572 
-0.658 


-0.181 
-0.086 


#7,C&N,t 
#7,ca,r 


ALS 11497,LSIV+10 9 


-4 50+ a62 

+' JU -0.57 


_-l n S +0.33 
J - -0.52 


-3.70 > 


0.06 


-0.241 


-0.158 


#7,ca,rv 


[SRA98] 6-1829 


-4.90 > 


-4.50 > 


-4.42 > 


0.03 


-0.830 


-0.074 


#7,ca 





PS 

r 

s 

S 



173153.70+064706.00 


sdB 


13 


.52 


13 


.74 


173630.26+073601.04 


sdB 


13 


.68 


13 


.91 


173651.20+280635.00A 


sdB 


1 1 


,53 


11 


,44 


173651.20+280635.00B 


F7V 










173812.47+263408.58 


sdB 


13 


.86 


14 


88 


175137.44+371952.37 


He-sdB 


13 


.25 


14 


,64 


175324.62+691058.08 


He-sdO 


13 


,84 


15 


46 


175340.57 500741.80A 


sdB 


13 


.21 


12 


88 


175340.57-50074 1.80B 


F7V 










175548.50+501210.77 


He-sdO 


13 


.31 


12 


.94 


181032.00+053909.00 


sdB 


13 


13 


13 


.81 


182530.60+084749.34 


sdO 


13 


.46 


14 


.28 


183129.61+085107.70 


sdO 


13 


.90 


14 


59 


184559.80 413827.00 


He-sdB 


13 


,89 


14 


.63 


190211.74 513005.61 


He-sdO 


11 


10 


11 


,53 


190302.50-352829.00 


sdB 


13 


,36 


14 


,34 


191109.30 140654.00 


He-sdO 


1 1 


.77 


11 


.77 


191109.30 140654.00 












191106.82+871604.35 


sdB 


12 


,86 


12 


.84 


200631.53+014307.25 


sdB 


13 


.74 


14 


45 


200951.66+031032.81 


sdB 


12 


,63 


12 


,63 


200918.25+095143.67 


sdB 


13 


,86 


14 


.58 


201554.27+054824.91 


sdB 


13 


,60 


14 


.41 


202027.26+070414.45A 


sdO 


13 


.85 


14 


,02 


202027.26+070414.45B 


F3V 










202059.85 225001.42A 


sdB 


1 1 


,92 


11 


.82 


202059. 85-225001.42B 


GOV 










202216.83+015225.18A 


sdB 


13 


,56 


13 


19 


202216.83+015225. 18B 


F6V 










202332.75+013618.82 


sdB 


12 


.78 


12 


.42 


203850.22-265747.16A 


sdO 


11 


.81 


11 


.90 


203850.22-265747.16B 


G3.5III 










203913.41+201309.38 


sdB 


12 


.69 


12 


.64 


204302.45+103427.08 


He-sdO 


1 1 


,62 


12 


,06 


204709.60+121015.38 


sdB 


13 


.73 


14 


.31 


205321.05 400626.94 


BHB 


11 


.55 


12 


.25 


205652.32+033548.88 


sdO 


13 


,36 


14 


.27 



27780+ 1 ™ 
28100^90 
36160+6|oo 
6240 

22880+ ™ 
34630+ soo 

43710+ 280 

+ J/1 -1160 

32430+ 88 °, 

6180 

40370! 
31890+ «o 

48970+ 800 

62840+f|}° 

35930+ 840 

72 3 00+ 5380 
'" u -3260 

32100+J™ 
55970+«40 
51300+«™ 

31610+l|g 

29600! 

27160+ 350 
4,iuu - 250 

27000! 340 
32510! 230 

442 10+ 2630 
++^ lu -2280 

6070 

34350! 14 o 
5920 

36700+ 11 60 
6310 

+ 460 



29450! 4 * 

58450! 4 ™ 
5270 



30600! 270 
45700! ' 6 °™ 
35510+ 300 
23050+ 
4035Q -n 3 °o 



5.43+ 



^ -25+0.18 
J - -0.07 
1+0.07 
-0.06 
c 94+0.84 
J - z +-0.84 

4.00 

5 56+ 05 

5 89+ 012 
J -° ? -o.io 

c. o,+0.41 
5 - S1 -0.06 
c qc+0.18 
-0.18 

4.04 

5 96+° 05 

J - ?u -0.15 
c qi+0.11 
5.81_ 010 

j+0.08 



5.89 



-0.11 



6.27+ 



,+0.19 
-0.08 
5 9^+0.27 

c 4 n+0.11 

;+0.31 
-0.30 

5 69+ 71 

-0.09 

5 50+ 33 
J - -0.00 



5.26+ 



5.81 + 



1+0.13 
-0.07 

5 62+ 06 

J - oz -0.06 

5 46+ 04 

J -+ u -0.05 

5 60+° 06 

J - uu -0.08 
c ss+005 
3-»»_0.06 

5 98+° 25 

J - 7o -0.23 

3.08 

S 44+019 
J -++-0.02 

4.04 

5 70+0.25 
J -' ? -0.12 

4.27 

5 60+° 06 

J - uu -0.04 

5 04+° 39 

J -"+-0.17 
3.02 

c 59+0.08 
J - J -0.06 
6.08+ 013 



-0.2 
c sc+0.07 

j+0.05 
-0.09 

5 55+ 07 
J - JJ -0.09 



4.78+ 



-2.53 > 
.9 10+0.34 

•'■"-0.78 

-1 09+° 69 



-3.26 > 

A 9O+0.04 

"°- 28 -0.06 

05+ 2 02 

u - UJ -0.02 
5+0.31 
-1.04 



-2.25+ 



1.30+ 



-v+0.27 
-0.10 

-2.46 > 



-2.52 



+0.29 
-0.29 



_9 4n +0.26 
i -+°-0.29 

2 10+' 10 
^- lu -0.38 
q 02+2.10 
-0.03 

-1.96 > 
25+ 085 



0.01 



-0.09 



-1.83+ 



2+0.09 
3 -0.09 
9 09+O.2I 
" Z - S2 -0.33 

-2 96+ 018 

■'■'"-0.26 

-2 54+°- 17 
i - J -0.20 

1 70+O.O6 
'■"-0.07 

_9 1 2+0.34 
i.ij_ 079 



-i-sotJo 



-1.52 

-2.44 
-1.13 



+0.22 
0.21 

+0.07 
0.20 

+0.27 
0.29 



-3.18 > 

0.78!°i 2 
, C9+0.05 

■ -0.09 

-1 36+° 09 

-0.09 
-2 98+ 024 



Table 2. Atmospheric parameters of hot subdwarf stars in the GALEX survey. Continued. 



GALEXJ 



Type 



Mjv V 
(mag) (mag) 



T e tr 
(K) 



(cm s -2 ) 



Abundances, [X/H], or [Fe/H] 



E(S -V) V - J J - H Comments Other name 



He 





"NT 

IN 


y J 


(mag) 


(mag) 


(mag) 






9 7O+0.20 
-3.61 > 


-2 90+ - 39 
-0.45 

-3.20 > 


-2.55 > 

-3.76 > 


0.02 
0.03 


-0.499 
0.611 


-0.145 
0.210 


#7,ca 
#7,r 


CD-43 14300 


[Fe/H] 


= 0.31 






(0.622) 


(0.173) 


M 




1/ty -0.18 

-2 46 +0 - 83 

^"-0.58 

-3 84+ - 28 
-0.36 


-2 05 +aM 

9 9 c +0.45 
-0.22 

-3 96+ 053 
-0.63 


-2 16+ 052 
-2.98 > 
-3.53 > 


0.03 
0.23 
0.06 


-0.711 
0.349 
-0.127 


-0.164 
0.220 
0.028 


#7,C&N 
#7,C&N 
#7,ca 




■3.97 > 


-4.19 +l ?-« 

—i.ZI 


-3.82 > 


0.07 


-0.511 


0.002 


s 


HS 2126+8320 


3.80 > 


-3.63 > 


-3.89 > 


0.01 


0.599 


0.233 


#7,r 




[Fe/H] 


= -0.48 






(0.552) 


(0.160) 


M 




-2 60+ 035 

z - DU -0.88 


-1.94 > 


-2.02 > 


0.08 


-0.296 


-0.161 


#7,C&N 




■3.79 > 


-3.25 > 


-3.36 > 


0.01 


0.557 


0.255 


#7,r 




[Fe/H] 


= -0.65 






(0.511) 


(0.244) 


M 




-j.Dl > 


—j. 04 > 


1 77 
— J. I 1 > 


U.Ul 


n A77 


n tit 


til r r\r 

tt 1 ,r,rv 




Ire/HJ 


— — U.JJ 






(U.JOJ ) 


(U. 1 l\J) 


JV1 




■5.07 > 


-4.20 > 


-4.41 > 


0.11 


-0.498 


-0.171 


#7,ca 


HS 2149+1428 


-3.06 > 


9 07+0.46 

J '-0.90 


-2.62 > 


0.09 


-0.027 


0.100 


#7 


PG 2158+082 


■4 62 +0 - 64 

-0.53 

-4.65 > 


A «+0.51 
^""-0.38 

-4.09 > 


-4 55+ 055 
-3.63 > 


0.02 
0.07 


-0.527 
0.277 


-0.036 
0.319 


#7 
#7 


BPS CS 30337-0074 
PG 2223+171 


■5.98 > 


-3.14 > 


-3.30 > 


0.00 


0.558 


0.100 


#7 


TYC 1703-394-1 


[Fe/H] 


= 0.08 






(0.701) 


(0.180) 


M 




-4.68 > 


_ 4 Q9+0.38 
■ -0.86 


-4.09 > 


0.02 


-0.593 


-0.109 


#7 


HS 2233+1418 


9 tc+0.45 
-0.19 


-2.45 > 


-2.03 > 


0.00 


-0.985 


-0.091 


#7,C&N 




-4.90 > 


-4.40 > 


-4.32 > 


0.00 


-0.489 


-0.181 


#7,ca 


Balloon 90900004 


-4.05 > 


—4 3Q+0-39 
H - J '-0.61 


-4.36 > 


0.04 


-0.214 


-0.169 


#7,s 




-3.98 > 


-3.86 > 


-3.71 > 


0.05 


-0.609 


-0.157 


na,f,t 


HS 2242+3206 


-4.65 > 


-4.51 > 


-4.54 > 


0.02 


-0.745 


-0.151 


#7,r,rv 




-4.88 > 


-4.54 > 


-4.33 > 


0.02 


-0.495 


-0.178 


na 




-4.48 > 


-4 46+ - 41 

-2.69 


-3.98 > 


0.03 


-0.544 


-0.187 


#7,s,rv 


PG 2321 + 142 


-2.64 > 


-2.53 > 


-3.08 > 


0.01 


0.850 


0.146 


na.he.r 




[Fe/H] 


= -0.29 






(0.685) 


(0.180) 


M 




-3.92 > 


-3.43 > 


-3.42 > 


0.00 


0.265 


0.179 


d 




[Fe/H] 


= -0.31 






(0.232) 


(0.119) 


M 




4.83 > 


-4 83+ - 37 


-4.23 > 


0.02 


-0.631 


-0.103 


#7,s 


CD-35 15910 


-3.97 > 


-4.08 > 


-3.99 > 


0.09 


-0.330 


-0.116 




FBS 2347+385 



205956.79-423203.58 


BHB 


1 1 


.54 


1 1 


,43 


210031.71+145213.83A 


sdB 


13 


.51 


13 


,09 


210031.71 + 145213.83B 


F5.5IV-V 










211041.43+142559.45 


He-sdO 


13 


16 


14 


.03 


211230.92 152410.37 


He-sdO 


1 1 


.46 


12 


,56 


212240.41+055808.08 


sdB 


13 


.87 


14 


.07 


212328.22+833351.21 


sdB 


13 


,84 


13 


,95 


212424.35+150619.44A 


sdB 


12 


.83 


13 


.06 


212424.35+150619.44B 


F7V 










212544.48+202812.18 


He-sdO 


12 


.21 


13 


67 


213730.94+221908.04A 


sdB 


13 


.03 


13 


.51 


213730.94+221908.04B 


G1V 










214022.87 371414.05A 


sdB 


12 


.48 


12 


.79 


214022.87-371414.05B 
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J1831+0851 is possibly an imminent WD progenitor. It lies 
well below the HeMS at a gravity higher than typical for He-sdO 
stars. Also, its He abundance is not comparable to its spectral class. 

J1845-4138. We found T eff = 35 930+ 4 **> K, log g = 5.23+°g 
and [He/H], = 2.1+H in good agreement with our previous analy- 
sis. This star shows an extremely high He abundance compared to 
its location in the r c ff-log g diagram. It is in the region of post-EHB 
evolution, where He-weak sdB stars are, but shows a He abundance 
typical for He-sdO stars. Its surface gravity is also inconsistent with 
its location. A high-dispersion spectrum is needed for further anal- 
ysis. No infrared excess was found. 

J1902-5130 (CD-51 11879, LSE 263) was analysed by Hus- 
feld et al. (1989). They found T eff = 70000 ± 2500 K, logg = 
4.9 ± 0.25, [H/He]. < -1.0, [C/He]. = -4.0 ± 0.5 and [N/He]. = 
-2.4 ± 0.3. Our analysis confirms the temperature and the He dom- 
inated, N-rich atmosphere, but provided a 0.5 dex higher surface 
gravity. Missing or incomplete atomic data and broadening param- 
eters in the earlier analysis can be accounted for such discrepancy. 
Due to upper limit determinations and asymmetric errors, a direct 
comparison with our results is not possible. 

J2043+1034 (LSIV+10 9) was analysed by Dreizler (1993) 
who found r eff = 44 500 ± 1000 K, logg = 5.55 ± 0.15 and mass 
fractions (j8 x = m x /m tBtd ): log/3 Hs . = -0.001, log/? c = -1.85 ± 0.1, 
log/? N = -2.45 ± 0.15 and log/So = -2.40 ± 0.15. Our analy- 
sis resulted in T cff = 45 700+ 1 ™ K, logg = 6.08^^, and mass 
fractions: log/? H e = -0.023, log/3 c = -2.69, log/3 N = -2.71 and 
logySo < -2.12. We derived a notably higher gravity and lower 
abundances, in particular for He and C. This can be due to the in- 
complete composition of our models and that we fitted the entire 
spectrum not just selected lines. The relatively lower quality of our 
fits for He-sdO stars also suggests the importance of heavier metals 
in their modelling. 

J2245+3221 (HS 2242+3206) was fitted by Edelmann et al. 
(2003). Their analysis of H r H 10 , Hen 4686, Hei 4471 and Hei 
4026 lines provided r cff = 29 300+ 1000 K, logg = 5.65±0.15 and 
[He/H], = -2. 8±0. 2. Our data covered the optical range from 3700 
to 7150 A and we derived T eff = 29 530^™ K, logg = 5.61^ and 
[He/H], = — 3 . l^o m g°°d agreement. 



9 CONCLUSIONS 

We presented a homogeneously modelled sample of 124 sdB and 
42 sdO stars and determined non-LTE atmospheric parameters by 
considering H, He and CNO opacities. With binary decomposition 
we were able to derive accurate atmospheric parameters for the 
members in ~30 composite spectra binaries. Our study shows that 
sdB stars concentrate in two groups in the T c g - log g and r cff - He 
diagrams, suggesting two typical H envelopes with different masses 
and compositions. He-sdO and He-weak sdO stars also clearly sep- 
arate in the r cff - He diagram and He- weak sdO stars can be related 
to sdB stars based on their He abundance and binary frequency. 

The developed method for spectral decomposition can reliably 
determine atmospheric parameters in composite spectra binaries. 
The method can also be used to remove the spectral contamination 
caused by merging stars from subdwarf spectra. Such cleaning pro- 
cess and homogeneous modelling on large samples are necessary 
to provide accurate atmospheric parameters for evolutional, pul- 
sational and diffusion theories. We investigated the impact of bi- 
narity on subdwarf atmospheric parameters and showed that unre- 
solved composite spectra binaries can cause significant shifts. Such 
systematic shifts can easily overwhelm non-LTE or abundance ef- 
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fects and we consider binary decomposition necessary for precision 
spectroscopy of hot subdwarfs with a notable IR excess. 

We did not find a significant sdB-A binary frequency with our 
current selection criteria. However, a joint analysis of ultraviolet 
and optical measurements would be necessary to spectroscopically 
distinguish the components in such binaries. Recently, Girven et 
al. (2012) investigated the companions of hot subdwarfs and found 
that F, G and K type MS companions outnumber A and M-type 
stars. 

Our abundance analysis of He-sdO stars confirms the C&N 
classification scheme of Stroeer et al. (2007): cooler stars show 
N enriched atmospheres, hotter He-sdO stars are C-rich while we 
detected both C and N mostly in between these two groups. The 
C&N-class objects suggest a C/N ratio of ~3 in He-sdO stars. The 
spectra of cooler sdB stars show a N enrichment while hotter sdB 
stars do not show notable CNO abundances in general. 

The observed luminosity distribution suggests a shorter EHB 
(~100 Myr) and longer post-EHB (~ 50 Myr) evolution than pre- 
vious studies found. The luminosity distribution of sdB stars is 
sharply peaked at log(L/L G ) » 1.5 assuming an average sdB mass 
of 0.48 M Q . We found two peaks in the luminosity distribution of 
sdO (post-EHB) and He-sdO stars. 

We reviewed our results in terms of the current evolutional 
theories. Our study does not give satisfying answers for the rela- 
tive contributions of the various channels, but shows evidences that 
He-sdO and He-rich sdB stars are related. The He-weak sdO se- 
quence and the group of He-rich sdB stars are mixtures of different 
populations. 

A binarity analysis of the hottest He-sdO stars (T eff > 60 000 
K) could provide interesting results about the relation of these stars 
to cooler He-sdO or He-weak sdO stars. In this region the se- 
quences of evolved He-sdO stars and He-weak sdO stars overlap. A 
much larger sample could help clarifying the surface gravity diver- 
sity observed at these stars, which currently suggests a connection 
to the He-weak sequence. 

Further follow-up is necessary at higher resolution and high 
signal-to-noise to improve CNO determination in the optical and 
better understand subdwarf stars, and the contributions of various 
formation and evolution channels. With higher quality data it may 
become possible to derive reliable O abundances that could offer 
further clues in understanding these stars. 

In addition to subdwarfs, our survey recovered six WDs that 
will be discussed in a forthcoming paper. 



10 THE SD1000 COLLABORATION 

A similar study on a larger sample could provide a more complete 
picture of hot subdwarfs. For this reason we would like to extend 
this pilot study and initiate a cooperative work under the name 
"SD1000 Collaboration", dedicated to process ~1000 hot subdwarf 
spectra in a homogeneous way. Such an analysis can provide im- 
portant observational constraints for theoretical studies and help the 
advance of subdwarf research. The Subdwarf Database currently 
holds 2399 stars and about 1900 of these have V, J and H photo- 
metric measurements. Figure 21 shows the colour indices of these 
stars. About ~800 stars in the data base are brighter than V = 15 
and ~ 1 150 are brighter than V = 16 providing a rich source of 
bright targets for the project. 

However, recording new, high-quality optical spectra (R > 
3000 or 8 A < 3 A, SNR > 100) for -1000 stars would require about 
60-80 nights on 4-m class telescopes what is not practically feasi- 
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Figure 21. V—J - J-H diagram for ~1900 hot subdwarf stars in The Subd- 
warf Database (0stensen 2006), data were collected from various sources. 
Typical error bars for colour indices were calculated from average photo- 
metric errors. 

ble. Acquiring data from public archives could be a way, although, 
data reduction of various observing runs can introduce a significant 
overhead and not all data are available in such archives. Moreover, 
it is important to maximise the number of spectra to number of in- 
struments ratio in the sample to decrease instrumental biases. The 
most efficient strategy would use contributed data for the bulk and 
obtain new observations for a small subset to fill the sample. There- 
fore, the success of this project depends strongly on contributed 
data and we look for collaborators. We estimate that enough ob- 
servations (~2000 stars) have been taken over the last 30 years to 
compile the target sample. 

Model atmosphere analysis would require ~2000 hours with 
single star models and 8 processors. Using binary decomposition 
and estimating about 20 per cent composite spectra in the sample 
would extend the analysis to ~2400 hours, which is still manage- 
able. 
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Figure Al. Panels a and d of Figure 9 complemented with error bars and upper limit measurements. Left: Observed C abundances with respect to effective 
temperature. The vertical line at 38 000 K separates sdB and sdO stars and the horizontal line shows the solar C abundance. Right: C and He abundance 
correlations. 
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Figure A2. Panels b and e of Figure 9 complemented with error bars and upper limit measurements. Left: Observed N abundances with respect to effective 
temperature. The vertical line at 38 000 K separates sdB and sdO stars and the horizontal line shows the solar N abundance. Right: N and He abundance 
correlations. 
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Figure A3. Panel c of Figure 9 complemented with error bars and upper limit measurements shows our observed O abundances with respect to effective 
temperature. The vertical line at 38 000 K separates sdB and sdO stars and the horizontal line shows the solar O abundance. 
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Figure A4. Panel /of Figure 9 complemented with error bars and upper limit measurements showing the C and N abundance correlations found in the GALEX 
sample. 



